/HO  2- 


AT> 


)?\ S' 


TECHNICAL  REPORT  NO.  11758 (LL141) 

AN  APPLICATION  OF  PLASTICITY  THEORY 
TO  THE  SOLUTION  OF  THE  RIGID  WHEEL- 
SOIL  INTERACTION  PROBLEM 

March  1973 


U.S.  ARMY 

TANK-AUTOMOTIVE 

COMMAND 

CONTRACT 

DAAE07-72-C-0033 


TECHNICAL  LIBRARY 


by 


LESLIE  L.  KARAFIATH* 
EDWARD  A.  NOWATZKI* 
I.  ROBERT  EHRLICH** 
JOHN  CAP  IN** 


GRUMMAN  AEROSPACE  CORPORATION* 
STEVENS  INSTITUTE  OF  TECHNOLOGY** 


MOBILITY  SYSTEMS  LABORATORY 


U.S.  ARMY  TANK  AUTOMOTIVE  COMMAND  Warren,  Mic 


Distribution  of  this  document  is  unlimited 


I 


Technical  Report  No. 11758  (LLl4l) 


AN  APPLICATION  OF  PLASTICITY  THEORY 
TO  THE  SOLUTION  OF  THE  RIGID  WHEEL  - 
SOIL  INTERACTION  PROBLEM 

MARCH  1973 


FINAL  REPORT 

by 

Leslie  L.  Karafiath* 
Edward  A.  Nowatzki* 
I.  Robert  Ehrlich** 
John  Cap in** 


Prepared  for  the 

United  States  Army  Tank -Automotive  Command 
Mobility  Systems  Laboratory 


Prepared  by 
Research  Department 
Grumman  Aerospace  Corporation 
Bet'hpage,  New  York  11714 

Under  Contract  No.  DAAE07-72-C-0033 


*  Grumman  Aerospace  Corporation 

**  Stevens  Institute  of  Technology 

Distribution  of  this  document  is  unlimited. 


FOREWORD 


One  of  the  most  urgent  needs  in  land  mobility  technology  is  for  a  valid 
comprehensive  methodology  to  (l)  evaluate  new  potentials  and  guide  their  exploi¬ 
tation;  (2)  support  complex  decision  processes  throughout  the  military  material 
development  cycle  with  analyses  which  incorporate  the  effect  of  land  mobility 
capabilities;  (3)  reduce  the  time  and  cost  required  to  develop  land  mobility 
systems,  in  particular  vehicles,  responsive  to  perceived  or  actual  needs. 

A  component  of  paramount  importance  of  the  required  land  mobility  technology 
is  modeling  or  simulating  the  interaction  of  the  terrain-vehicle -man  system. 

In  1971,  the  U.S.  Army  Tank -Automotive  Command  (TACOM)  and  the  U.S.  Army 
Corps  of  Engineers,  Waterways  Experiment  Station  (WES)  completed  a  first  gener¬ 
ation  terrain-vehicle-man  interaction  simulation  called  the  "AMC  *71  Vehicle 
Mobility  Model." 

Further  research  effort  is  needed,  however,  to  improve  the  accuracy  and 
range  of  applicability  of  this  model.  One  of  the  sub-models  to  be  perfected  is 
the  analysis  and  simulation  of  a  wheeled  vehicle  moving  in  soft  soil. 

With  this  in  mind,  TACOM  has  contracted  with  Grumman  Aerospace  Corporation 
to  develop  a  rigorous  computerized  scheme  for  calculating  rigid  wheel  slip, 
sinkage,  torque  requirements,  motion  resistance  and  drawbar  pull  from  vehicle 
and  soil  inputs,  such  as,  wheel  load,  width,  diameter,  soil  cohesion,  friction 
and  bulk  density. 

This  was  conceived  to  be  the  first  step  (the  second  being  a  similar  scheme 
for  pneumatic  tires)  toward  the  establishment  of  a  rigorous  method  for  predicting 
the  maximum  speed  that  a  wheeled  vehicle  can  attain  in  soft  soil  as  well  as  for 
the  accurate  assessment  of  performance  of  fuel  consumption  requirements  in  per¬ 
forming  specific  missions. 
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ABSTRACT 


Plasticity  theory  is  applied  to  the  analysis  of  soil-wheel  interaction  . 

The  problem  is  reduced  to  the  determination  of  stresses  at  the  interface  of  the 
rigid  wheel  and  the  soil.  Once  these  stresses  are  known,  the  wheel  load,  torque, 
pull  and  drag  are  obtained  by  integrating  the  stresses  along  the  wheel  perimeter. 

To  find  the  distribution  of  interface  stresses,  the  basic  differential  equations 
of  equilibrium  are  combined  with  the  Coulomb  -  Mohr  yield  criterion  for  soils  and 
the  equations  are  solved  by  a  numerical  procedure.  The  numerical  solution  scheme 
and  the  computer  program  which  accomplishes  the  solution  are  described  in  detail. 

Tests  performed  with  the  soil  bin  dynamometer  facilities  of  Stevens  Insti¬ 
tute  of  Technology  are  discussed.  Test  results  show  good  agreement  with  prediction. 
The  validity  of  assumptions  introduced  in  the  computational  scheme  is  examined. 

It  is  found  that  refinement  in  the  assumptions  regarding  the  distribution  of  inter¬ 
face  friction  and  the  magnitude  of  the  "separation  angle"  would  further  improve 
the  accuracy  of  the  method. 

Finally,  it  is  concluded  that  the  proposed  theory  for  predicting  rigid 
wheel  performance  is  fundamentally  correct  and  is  practical  from  the  viewpoint  of 
required  computer  time. 

This  report  represents  an  essential  first  step  toward  the  establishment  of 
a  rigorous  simulation  of  the  soft  soil  performance  of  wheeled  vehicles. 
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NOTATION 
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=  singular  point 
=  cohesion 
=  constants 

=  drawbar  pull,  drag 
=  function 
=  slip 

=  constant  defining  threshold  slip 

=  constant  in  plate  sinkage  equation 
=  constant  in  slip  equation 
=  load 

=  exponent  in  plate  sinkage  equation 
=  bearing  capacity  factor 
=  pressure  in  plate  sinkage  equation 
=  constant  in  plate  sinkage  equation 

=  normal  stress 

=  normal  stress  in  forward  field  at  a 

m 

=  normal  stress  in  rear  field  at  ot 

m 

=  normal  stress  at  the  edge  of  the  wheel  corresponding  to  transverse 
failure 

=  radius  of  wheel 
=  (cr^  +  a  )/ 2,  reference  stress 

=  geometric  coordinates 
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a 


=  central  angle  (measured  from  vertical) 

o'  =  entry  angle 

0 

<¥  =  angle  of  separation 

=  rear  angle 

Y  =  unit  weight  of  soil 

6  =  angle  of  inclination  of  resultant  stress  to  normal,  angle  of 

shear  mobilization 

6  =  interface  friction  factor  for  cohesive  soils 

c 

e  =  slope  angle 

9  =  angle  between  x  axis  and  major  principal  stress 

X  =  load  coefficient  (l/Lq) 

h  =  tt/4  -  cp/2 

§  =  tolerance  limit 

o  =  (ct1  +  ct3)/2  +  ♦ 

cr  =  normal  stress 

n 

a  =  principal  stresses 

1,  o 

t  =  shear  stress/shear  strength 

t  =  maximum  available  shear  strength 

max 

Tmob  =  s-iear  strength 

cp  =  angle  of  internal  friction 

ik  =  c  cot  cp 

(jo  =  pull  coefficient*  (DB/L  ) 

SUBSCRIPTS 

a  =  act ive 

i,  3  =  designation  for  slip  line  numerals 

K  =  designation  for  iteration  numerals 

m  =  model 

max  =  maximum 

min  =  minimum 

mob  =  mobilized 

=  input  (design)  values 
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I.  SCOPE  OF  WORK 

The  scope  of  work  as  described  in  the  Request  for  (Contract)  Proposal  work 
statement  was  to  conduct  theoretical  and  experimental  investigations  on  the 
application  of  plasticity  theory  to  the  solution  of  wheel  soil  interaction  pro¬ 
blems.  Emphasis  was  placed  on  the  following  items: 

•  Development  of  a  computer  program  for  the  prediction  of  performance 
of  rigid  wheels 

•  Performance  of  validation  tests  to  compare  theoretical  and  experimental 
results 

Presentation  of  the  results  of  the  research  work  is  arranged  in  the 
above  order  with  an  introductory  section  on  theoretical  background. 
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II.  THEORETICAL  BACKGROUND 


A.  Basic  Assumptions  and  Definitions 

In  the  application  of  plasticity  theory  to  the  solution  of  the  soil- 
wheel  interaction  problem,  the  following  basic  assumptions  are  made: 

•  The  soil  is  semiinfinite,  homogeneous,  and  isotropic. 

•  Strength  properties  of  the  soil  are  characterized  by  the  cohesion 
and  internal  angle  of  friction. 

•  Soil  failure  is  governed  by  the  Mohr-Coulomb  yield  (failure)  criterion. 

•  Stresses  are  the  same  in  all  planes  parallel  to  the  plane  of  notion 
of  the  wheel  (the  problem  is  assumed  to  be  two  dimensional). 

•  Wheel  velocity  is  constant. 

•  Soil  inertia  forces  are  negligible. 

•  Pore  water  pressures  are  negligible. 

Some  of  the  terms  frequently  used  in  the  following  discussions  are 
defined  below: 

•  Plastic  state  of  stress,  plastic  equilibrium  conditions,  or  failure 
conditions  are  used  interchangeably  to  define  a  stress  state  at  any 
point  in  the  soil  where  the  Mohr-Coulomb  yield  criterion  is  satisfied. 

•  Slip  line  is  a  line  along  which  failure  conditions  obtain. 

•  Slip  zone,  slip  line  field,  or  failure  zone  is  a  finite  area  where 
the  Mohr-Coulomb  yield  criterion  is  satisfied. 

•  Angle  of  interface  friction  is  the  angle  that  determines  the  shear 
stress  in  a  Mohr-Coulomb  plot  (Fig.  l)  according  to  the  following 
relationship 

T  =  ^an  +  c,co^cP)^an  ^  (l) 

If  6  =  cp,  Eq.  (l)  reverts  to  the  Mohr-Coulomb  yield  criterion. 
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B.  Soil-Wheel  Interaction  Concept 

The  application  of  plasticity  theory  to  the  analysis  of  soil-wheel 
interaction  is  based  on  the  concept  that  failure  conditions  develop  in  the 
soil  beneath  towed  or  driven  wheels,  and  that  they  control  the  interface 
stresses.  The  validity  of  this  concept  is  supported  by  experimental  obser¬ 
vations  as  well  as  by  the  good  agreement  between  measured  interface  stresses 
and  those  predicted  by  the  theory.  The  concept  and  its  practical  application 
for  wheel  performance  calculations  is  dealt  with  in  detail  in  References  1 
and  2  and  is  discussed  here  only  briefly  for  the  convenience  of  reference. 

Any  soil-wheel  interaction  problem  can  be  reduced  to  the  determination 
of  interface  stresses.  Once  these  are  known,  the  wheel  can  be  considered  as 
a  free  body  and  the  computation  of  load,  torque,  and  drag  or  drawbar  pull 
consists  of  a  simple  integration  of  normal  and  shear  stresses  along  the  wheel 
perimeter. 

In  soft  soils,  the  interface  stresses  are  controlled  by  the  inability 
of  the  soil  to  carry  loads  higher  than  those  causing  failure  in  the  soil. 
Failure  in  the  soil  occurs  in  zones  along  slip  lines.  There  are  generally 
two  failure  zones  beneath  loaded  wheels,  as  shown  in  Fig.  2,  a  forward  and 
a  rear  one.  Each  failure  zone  is  comprised  of  three  different  types  of 
failure  states,  as  described  in  Fig.  2.  The  adjectives  denote  both  the 
location  of  these  zones  and  the  direction  in  which  failure  occurs.  One 
requirement  is  that  the  two  zones  must  meet  and  that  at  their  common  point 
the  interface  stresses  be  the  same.  In  certain  cases,  there  is  only  a  single 
failure  zone,  as  shown  in  Fig.  3-  This  occurs  when  the  interface  stress  for 
the  single  zone  is  less  than  that  which  would  cause  failure  in  the  other 
direction  as,  for  example,  at  point  B  in  Fig.  3. 

C.  Interface  Stresses  and  Wheel-Soil  Interaction 

The  role  of  the  wheel  in  soil-wheel  interaction  is  twofold.  First,  the 
entry  and  rear  angles  are  the  limits  of  the  integration  of  interface  stresses 
that  yield  the  load,  torque,  and  drawbar  pull.  The  integration  limits  express 
the  change  in  the  contact  area  and  thus  affect  wheel  performance  profoundly. 
Second,  the  applied  torque  on  the  wheel  generates  shear  stresses  at  the 
interface  that  influence  the  magnitude  and  distribution  of  normal  stresses 
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that  the  soil  can  carry.  The  entry  and  exit  angles  also  affect  the  failure 
stresses  in  the  soil.  These  interaction  effects  are  illustrated  in  Figs.  4 
and  5»  The  method  of  calculation  of  slip  line  field  geometries  and  inter¬ 
face  stresses  is  described  in  Section  III;  here  the  results  are  anticipated 
for  the  purpose  of  realistic  illustration. 

Rear  slip  line  fields  for  various  angles  of  interface  friction  angle 
are  shown  in  Fig.  4.  The  rear  field  was  selected  for  illustration  because 
of  the  drastic  effect  of  the  interface  friction  angle  on  the  geometry  of  the 
rear  slip  line  field.  For  comparison,  all  parameters  but  the  angle  of  inter¬ 
face  friction  are  the  same  in  Fig.  4.  The  diminishing  of  the  extent  of  the 
rear  slip  line  field  with  the  increase  of  the  interface  friction  angle  is  evident 
from  Fig.  4.  With  the  change  in  the  geometry,  the  associated  normal  interface 
stresses  also  change  as  shown  on  the  right  side  of  Fig.  4.  The  effect  of  the 
angle  of  interface  friction  (generated  by  the  applied  torque)  on  the  normal 
stresses  is  very  pronounced  as  seen  in  the  illustration. 

The  effect  of  the  entry  and  exit  angles  on  the  interface  stresses  is 
less  pronounced.  Slip  line  field  geometries  for  various  entry  angles  are  com¬ 
pared  in  Fig.  5j  all  other  variables  being  the  same.  The  change  in  the  extent 
of  slip  line  fields  is  due  mostly  to  the  condition  that  they  all  were  made  to 
end  at  the  angle  of  15  degrees.  The  corresponding  normal  stresses  shown  on  the 
right  side  of  the  illustration  appear  to  rise  along  parallel  lines;  however, 
there  is  a  significant  difference  in  the  rate  of  the  rise  as  a  comparison  of  the 
dashed  line  with  the  solid  one  at  60°  shows.  The  dashed  line  is  the  normal 
stress  for  30°-entry  angle  but  plotted  at  60°  so  that  a  comparison  with  the 
solid  line  for  60°-entry  angle  be  facilitated.  This  comparison  shows  that  the 
entry  angle  influences  the  normal  stresses  but  to  a  lesser  degree  than  the 
angle  of  interface  friction. 

D.  Application  of  Plasticity  Theory  to  the  Determination  of  the  Geometry  of 

Soil  Failure  Zones  and  Associated  Interface  Stresses 

In  plasticity  theory,  the  differential  equations  of  equilibrium  are  com¬ 
bined  with  a  yield  criterion  for  the  material.  For  soils,  the  failure  criterion 
that  is  generally  recognized  to  give  the  best  agreement  with  experiments  is  the 
Mohr-Coulomb  yield  criterion,  which  expresses  the  following  linear  relationship 
between  yield  shear  strength  and  normal  stress  (for  notations  see  page  vi): 


(2) 


t  ss  c  +  o  tan  cp 
n 


This  failure  criterion,  when  combined  with  the  differential  equations 
of  equilibrium,  yields,  after  some  manipulations  (Refs.  3  -  5),  the  following 
set  of  differential  equations: 

dz  =  dx  tan  (9  ±  p,) 


V 

dn  ±  2o  tan  cpd9  =  -  sin  (e  ±  cp)  dx  +  cos  (e  ±  cp)  dz 

cos  cp  L  T  T  _ 


(3) 


The  cohesion  term  (c)  in  Eq.  (2)  does  not  appear  explicitly  in  Eqs.  (3); 
it  is  included,  however,  in  the  value  of  cr  which  equals 

c  =  l/2  (ct^  +  CT3)  +  c  cot  cp  (4) 

Equations  (3)  lose  their  meaning  if  cp  =  0,  since  cr  becomes  infinite. 

In  this  case,  the  following  differential  equations  apply  (Ref.  3): 

dz  =  dx  tan  (9  ±  p) 

ds  ±  2sd9  =  y  dz  (5) 

where  s  =  l/2  (o^  +  ^) 

For  properly  defined  boundary  conditions,  the  above  sets  of  the  differ¬ 
ential  equations  of  plasticity  for  soils  yield  a  unique  solution  in  the  form 
of  a  slip  line  field  and  associated  stresses.  At  any  point  within  the  slip 
line  field,  the  applicable  differential  equations  set  forth  above  are  satisfied. 

There  are  some  fundamental  aspects  of  the  theory  that  need  to  be  empha¬ 
sized  here  before  proceeding  to  the  discussion  of  the  numerical  procedures 
applied  for  the  solution  of  the  above  differential  equations.  These  are: 


•  The  Mohr-Coulomb  yield  criterion  as  applied  in  the  above  equation 
refers  to  effective  stresses,  i.e.,  for  the  case  where  pore  water 
pressures  are  negligible.  (Apparent  cohesion  due  to  pore  water  tension 
may  be  considered  as  effective  stress.  ) 

•  Equations  (3)  of  the  plasticity  theory  are  valid,  however,  for  a  non¬ 
linear  strength  envelope  (Ref.  6)  and  may  be  expanded  to  include  pore 
water  pressures  (Ref.  2)  or  soil  inertia  forces  (Ref.  7)» 
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•  The  Mohr-Coulomb  yield  criterion  implies  that  the  soil  strength  ex¬ 
pressed  hy  Eq.  (2)  is  available  regardless  of  the  volumetric  strain 
that  is  associated  with  the  stress  state  expressed  by  Eq.  (2).  In 
soil-wheel  interaction  problems  in  soft  soils,  where  the  soil  is  pro¬ 
gressively  compressed  as  the  wheel  advances,  the  volumetric  strain  in 
the  soil  is  generally  much  larger  than  in  static  problems  and  the 
shear  strength,  as  expressed  by  Eq.  (2),  is  available  for  a  failure 
zone  to  form.  However,  at,  and  in  the  close  vicinity  of  the  entry 
angle,  the  volumetric  strain  may  not  be  enough  to  mobilize  the  full 
Mohr-Coulomb  strength  and,  consequently,  the  rise  of  normal  stresses 
in  cohesive  soils,  though  rapid,  is  not  instantaneous  as  predicted  by 
the  theory. 

•  In  plasticity  theory,  solutions  obtained  by  integration  of  Eqs.  (3) 
are  termed  statically  admissible  solutions  and  are  considered  lower 
bound  solutions.  Kinematic  admissibility  is  analyzed  by  constructing 
velocity  fields  for  the  slip  zones  on  the  basis  that  the  material  is 
incompressible.  A  kinematically  admissible  solution  would  constitute 
an  upper  bound.  For  soil-wheel  interaction  problems  in  soft  soils, 
the  assumption  of  incompressibility  is  inappropriate  and  conclusions 
drawn  on  the  basis  of  such  an  assumption  are  inapplicable.  Experimental 
evidence,  as  discussed  in  Section  IV,  shows  that  interface  stresses 
predicted  by  the  application  of  Eqs.  (3)  agree  well  with  measured  ones 
indicating  the  validity  of  the  solutions. 

E.  Similitude  and  the  Application  of  Plasticity  Theory  to  Soil-Wheel 

Interaction 

The  concept  that  interface  stresses  are  governed  by  the  differential 
equations  of  plasticity  for  soils  leads  to  interesting  considerations  regard¬ 
ing  the  use  of  the  principles  of  similitude  in  soil-wheel  interaction  studies. 
Equations  (3)  can  be  written  in  dimensionless  form  if  a  characteristic  length 
and  a  reference  stress  is  introduced.  For  soil-wheel  interaction  problems,  the 
characteristic  length  may  be  conveniently  chosen  as  the  radius  of  the  wheel. 

The  geometry  of  the  slip  line  fields  in  a  soil-wheel  interaction  problem  will 
be  similar  for  the  same  cp  and  8  if  the  following  equality  holds:  . 
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(6) 
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The  subscripts  m  and  P  denote  model  and  prototype,  respectively. 
For  geometrically  similar  slip  line  fields 


<7 
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Thus,  at  the  interface,  the  distribution  of  c  values  are  the  same  but 
their  magnitude  changes  according  to  relationship  (7).  At  this  point, 
attention  is  called  to  the  definition  of  a  in  plasticity  theory  j^Eq.  (4)  . 

If  c  =  0,  or  cd  =  0,  and  Eqs.  (5)  are  used,  the  normal  and  shear  stresses  at 
the  interface  are  proportional  to  a.  As  a  consequence,  dimensionless  wheel 
performance  parameters,  such  as  pull  coefficients  (DB/L)  will  be  the  same 
for  all  geometrically  similar  slip  line  fields,  since  a  constant  ratio  of 
interface  stresses  can  be  factored  out  in  the  integration  for  both  load  and 
drawbar  pull.  Thus,  in  the  case  of  c  =  0  or  cp  =  0  soils,  the  application  of 
plasticity  theory  for  soil-wheel  interaction  confirms  the  validity  of  the 
results  obtained  in  similitude  studies  regarding  the  selection  of  dimensionless 
parameters,  at  least  for  wheels  traveling  at  velocities  where  soil  inertia 
effects  are  negligible. 

In  the  case  of  c  ^  0,  co  ^  0  (c  -  cp  soils)  geometrically  similar  slip 
line  fields  result  in  similar  distributions  of  a  along  the  interface.  How¬ 
ever,  in  this  case,  distribution  of  normal  stresses  is  no  longer  similar  to  the 
ct  distribution,  since  c  contains  the  term  c  cot  cp  (Eq.  4)  and  does  not. 

Thus,  similarity  exists  between  the  distribution  of  +  c  cot  cp,  the  shear 
g~fcp0gggs  t,  and  o.  In  the  integration  formulas  for  load,  drawbar  pull,  and 
torque,  the  constant  c  cot  cd  can  be  computed  separately  and  a  correction 
factor  for  c  -  cp  soils  in  the  similitude  relations  may  be  established. 

It  is  also  noted  that  while  Eqs.  (3)  also  apply  for  nonlinear  yield  cri¬ 
teria  (Ref.  6),  they  cannot  be  put  in  dimensionless  form  unless  the  yield 
criterion  is  linear.  Thus,  similitude  in  soil-wheel  interaction  as  expressed 
by  Eqs  (6)  and  (7)  is  restricted  to  soils  that  follow  the  linear  Mohr-Coulomb 

yield  criterion. 
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F.  Slip  and  the  Development  of  Interface  Friction 

For  slip  line  fields  to  be  determined  uniquely  by  the  differential 
equations  of  plasticity,  it  is  necessary  to  set  some  of  the  boundary  conditions 
at  the  interface.  In  this  so-called  mixed  boundary  value  problem,  two  variables 
of  the  four  independent  ones  in  Eqs.  (3)  must  be  specified.  In  the  solution 
for  the  conventional  bearing  capacity  problem,  z  and  9  are  specified  at  the 
base  of  the  bearing  plate.  In  the  soil-wheel  interaction  problem,  it  is  not 
possible  to  specify  the  value  a  priori  of  any  of  two  of  the  four  variables  at 
the  interface.  Instead,  the  boundary  conditions  take  the  form  of  a  relation¬ 
ship  between  x  and  z  given  by  the  wheel  geometry  and  a  relationship  between 
9  and  the  angle  of  interface  friction.  This  relationship  is  (Ref.  3) 


6  »  \  +  |  (A+6)  -  or 


(8) 


where 


A  = 


arc  sin 


Thus,  the  interface  friction  angle  needs  to  be  specified  along  the  inter¬ 
face  for  the  slip  line  field  and  the  associated  interface  stresses  to  be 
uniquely  defined.  In  wheel  performance  calculations,  however,  the  performance 
parameters  are  related  to  slip  rather  than  to  the  interface  friction  angle. 

The  development  of  shear  stresses  at  the  interface  is  associated  with  slip, 
and  mathematical  formulations  for  the  relationship  between  shear  stress  and 
slip  have  been  proposed  by  various  researchers.  On  the  basis  of  direct  shear 
tests,  Janosi  and  Hanamoto  proposed  the  following  relationship  between  mobilized 
shear  and  slip  for  tracked  vehicles  (Ref.  13): 


t  =  t 
mob  max 


(l  -  e-J/K). 


(10) 


For  compressible  soils,  which  are  of  primary  interest  in  off-road  loco¬ 
motion,  this  equation  properly  describes  the  relationship  between  shear  stress 
and  slip.  When  this  relationship  is  applied  to  the  rigid  wheel,  however,  it 
is  useful  to  include  a  constant,  in  the  slip  term  to  account  for  the  fact 
that  a  threshold  perimeter  shear  exists,  at  which  movement  of  the  wheel  starts. 
Thus,  Eq.  (10)  is  modified  as  follows: 

Vb  =  w((1  -  e'(3  *  3°,/K)  (11) 
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The  following  relationship  holds  between  the  shear  strength  mobilized  at  the 
interface  and  the  angle  6  (Fig.  l): 


tan  6  = 


mob 

a  +  ^ 

n 


(12) 


From  Eqs.  (ll)  and  (12)  comes  the  following  relationship: 


tan  6  =  tan  6  (1  -  e 

max 


( 1  -  e-(3+Jo>/K) 


(13) 


The  relationship  between  the  interface  friction  angle  and  slip  established 
by  Eq.  (13)  allows  the  computation  of  slip  for  various  values  of  6  if  j  and  K 
are  known.  The  concept  of  soil -wheel  interaction,  as  outlined  in  the  preceding 
paragraphs,  has  important  implications  regarding  the  value  of  6  .  According  to 

this  concept,  the  soil  adjoining  the  interface  is  in  the  active  state  of  failure. 
For  a  given  normal  and  shear  stress  at  the  interface,  there  is  one  Mohr  circle 
that  represents  the  active  and  another  one  that  represents  the  passive  state  of 
stresses.  Figure  6  shows  Mohr  circles  for  the  active  and  passive  state  for 
the  same  normal  stress  but  with  increasing  interface  shear  stress.  Stress 
circles  for  the  active  state  are  shown  by  full  lines,  and  by  dashed  lines  for 
the  passive  state.  The  interface  shear  stress  t  is  shown  to  increase  with  the 
interface  friction  angle  (t  >  T2  >  ti  and  §3  >  §2  In  the  active 

state,  the  center  of  Mohr  circles  is  to  the  left  of  the  shear  stress  ordinate; 
in  the  passive  state,  it  is  to  the  right.  From  the  construction  of  the  Mohr 
circles,  it  is  obvious  that  the  maximum  shear  stress  that  can  be  mobilized  in 
the  active  state  is  the  one  corresponding  to  a  Mohr  circle  that  has  its  center 
at  cr^  (circle  3  in  Fig.  6).  Were  the  shear  stress  higher  than  this,  the 
corresponding  Mohr  circle  would  represent  a  passive  state.  Thus,  it  is  in¬ 
correct  to  assume  that  the  full  soil  strength  can  be  mobilized  beneath  a 
wheel  or  track.  The  passive  state  beneath  a  wheel  or  track  can  exist  only  if 
the  soil  is  pushed  toward  the  wheel  or  track,  an  obviously  meaningless  situa- 
tion  for  vehicle  mobility.  The  only  possible  stress  state  in  the  soil  beneath 
a  wheel  or  track  is  the  active  state  of  stresses  and  it  follows  from  this  state 
that  the  maximum  mobilized  shear  stress  cannot  exceed  that  defined  by  Mohr 
circle  3  in  Fig.  6.  In  mathematical  terms 
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t  =  (<J  +  i|r )  tan  6 

max  n  max 


(CTn  +  If)  sin  cp 


(14) 


6  =  arctan  (sin  cp)  (15) 

max 


G.  Constraints  Imposed  on  the  Slip  Line  Fields  by  the  Wheel 

The  interface  friction  angle  degines  0  at  the  soil-wheel  interface  (Eq,  (8)) 
and,  together  with  the  geometry,  defines  the  boundary  conditions  so  that 
Eqs.  (3)  yield  a  unique  solution  for  a  particular  slip  line  field.  For  the 
soil-wheel  interaction  problem,  however,  the  boundary  conditions  for  both 
the  forward  and  rear  slip  line  field  must  be  defined  so  that  the  solution  be  uni¬ 
que.  This'  involves  the  determination  of  angle  or  that  defines  the  end  point  of  the 
forward  and  rear  field.  The  following  considerations  apply  in  this  respect. 

For  dry  sand  it  was  found  by  Sela  (Ref.  8)  that  the  angle  of  separation 
approximately  equals  the  developed  friction  angle.  This  finding  is  consistent 
with  the  concept  that  the  forward  failure  zone  extends  over  that  part  of  the 
wheel  perimeter  where  the  component  of  the  normal  and  shear  stresses  (AD)  in 
the  direction  of  motion  is  negative  (i.e.,  resisting  the  motion),  and  that 
the  backward  zone  extends  over  that  part  of  the  wheel  perimeter  where  this 
component  is  positive.  Applying  this  concept  to  soils  with  cohesion  and 
pore  water  pressures  results  in  the  following  relations: 


AD 

=  T  - 

cos 

a  - 

a  sin 

<y 

moo 

n 

Tmob 

=  (an 

+  \|r) 

tan 

6  (Eq. 

(12)) 

AD 

II 

+  if) 

tan 

6  tan  a  =  0 

/ 

fo  +  ^ 

\ 

<y 

=  arc 

tan  ( 

n 

—  tan 

§1 

m 


(16) 
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III.  COMPUTATIONAL  METHODS 


A.  General 

The  solution  of  differential  Equations  (3)  yields  the  slip  line  fields 
and  associated  interface  stresses  for  wheel  performance  calculations.  For 
numerical  computations,  differential  Equations  (3)  are  replaced  by  the 
following  finite  difference  Equations  (Fig.  7): 

i,0  i-l,0  i,I-l  1  i,0-l  2  i-l,o '  v  1  2' 

z.  .  =  z.  t  .  +  o',  (x.  .  -  x.  .  .)  (17 

1,0  1-1,0  2  1,0  1-1,0 


ct.  .  = 
i,I 


y(Cct.  +  Du  .)  +  2o.  o  .  fl  +  (9.  .  -  0.  .)  tan  cp 

_ -1-? j~i _ 1~1? 3 _ j~i  i~i? j  L _ j~i _ -l~i>j 

a  +  a 

i, 0-1  i-1,0 


o'.  .  -  o.  +  2  tan  cp(o.  0.  .  +  o  .  9.  .)  +  Y  (D-C) 

1,0  -1  1-1,0  1,0-1  1,0-1  1-1,0  1-1,0 


2  tan  cp  (o.  .  +  a.  . ) 

i,0-l  i-l,o 


where  x.  . ,  z.  .  are  coordinates  of  the  subscripted  nodal  point 
3-  ?  3  i  ?  3 

“l  ‘  tan  (0i,o-l  +  •‘t  °2  *  tan  (ei-x,3  -  >*> 


0  =  sln  (e-T)  (x.  .  -  X.  ,  .)  +  SSlM  (z.  .  .  «  .) 

cos  cp  1,0  1-1,0  COS  cp  1,0  1-1,0 


d  =  siL iiai  (x.  . .  x  .)  +  oosj^o}  (  .  _ . 

cos  cp  1,0  1-1,0  cos  cp  1,0  i-l,0 


These  difference  equations  permit  the  computation  of  the  coordinates  of 
a  nodal  point  (intersection  of  slip  lines),  as  well  as  the  values  of  o  and  0 
at  that  point,  from  the  known  values  at  neighboring  nodal  points  having  lesser 
subscripts. 
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The  slip  line  field  in  the  passive  zone  can  he  computed  by  equations 
starting  with  the  boundary  values  given  at  the  soil  surface  (sloping  or 
level).  In  the  radial  shear  zone,  the  same  equations  are  used,  but  special 
consideration  is  given  to  the  central  point  ("A")  where  the  second  family 
of  slip  lines  converge.  This  point  is  a  degenerated  slip  line,  where  © 
changes  from  the  value  at  the  passive  zone  boundary  to  that  specified  at  the 
active  zone  boundary.  The  total  change  in  0  is  divided  by  the  number  of 
slip  lines  converging  at  this  point  to  result  in  an  equal  A0  increment  be¬ 
tween  two  adjacent  slip  lines.  The  o  values  for  each  increment  are  computed 

2  (  Q  mm  O  \  *4-  p  VJ  fQ 

from  the  equation  a  =  e  K  o'  which  is  the  solution  of  the 

differential  equations  of  Eqs.  (3)  if  both  dx  and  dz  vanish.  With  these  © 
and  o  values  assigned  to  each  slip  line  at  the  singular  point,  the  coordinates 
as  well  as  the  o  and  0  values  for  all  other  points  in  the  radial  shear  zone 
can  be  computed  by  Eqs.  (17). 

In  the  active  zone,  the  same  equations  are  used  except  for  the  boundary 
at  the  interface,  where  ©.  .is  specified  and  x  and  z  must  lie  on  the  circle 

1 , 0 

with  radius  R.  For  numerical  computations,  the  circle  is  approximated  by  a 
polygon,  allowing  the  use  of  the  following  difference  equations: 


x.  .  = 


x.  .  .  +  o'  _F  x. 


i,j  1  +  Q'qF  i-1,0  0 


X  •  «,  .  _  +  Qf _  ( Z .  n  -  Z.  -  .) 

l-l, 3-1  0  l-l, o-l  i-l,o 


z.  . 

1  _ 


.  =  z.  .  .  _  +  F  (x.  _  .  ,  -  x.  .) 

,0  l-l, o-l  l-l, o-l  1,0 

o.  .  =  CT  .  +  o  tan  9  CT-  n  ■  +  (0.  .  -  0.  ,  .)  +  YC 

1,0  1-1,0  2  1-1,0  1,0  1-1,0 


where 


Q!q  =  cotan 


(2  (6i,o-l  +  0i,O}  “  ^ 


(18) 


F  =  tan  a. 

l-l, 0-1 


To  improve  the  accuracy  of  calculations,  an  iteration  is  performed 

where  6.  n  .,®.  .,,cr  .,a.  values  in  the  above  difference  equations 

i-l,3  i,o-l  1-1,0*  i,o-l 

are  replaced  by  values  averaged  with  the  computed  value  of  0.  .  and  o.  . , 

1,0  1,0 

respectively. 
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B.  Computation  of  a  Single  Slip  Line  Field 

The  finite  difference  equations  given  in  Section  III.  A  are  suitable  for 
the  numerical  computation  of  a  single  slip  line  field.  For  the  problem  of 
soil-wheel  interaction,  it  is  convenient  to  carry  out  the  computations  in 
such  a  way  that  the  slip  line  comprising  the  outer  boundary  of  the  field  ends 
at  a  specified  location  at  the  interface.  To  this  end,  a  sequence  of  opera¬ 
tions,  different  from  that  used  for  the  conventional  bearing  capacity  problem, 
is  employed,  as  shown  in  the  flow  diagram  in  Fig.  8.  Instead  of  computing  the 
variables  first  in  the  passive,  then  in  the  radial,  and  finally  in  the  active 
zone,  the  variables  are  computed  along  the  first  "j"  line  in  all  three  zones 
and  then  along  subsequent  "j"  lines  until  a  "j"  line  exceeds  the  end  point. 

Then  a  ",i"  line  is'  interpolated  so  that  it  ends  up  at  a  within  the 

m 

limits  of  tolerance  (?).  This  method  eliminates  the  time  consuming  trial  and 
error  procedure  of  finding  the  length  of  passive  zone  that  matches  the  arc 
length  of  the  active  zone.  The  length  of  the  passive  zone  in  this  procedure 
is  overestimated  so  that  the  last  "j"  line  overshoots  the  or  angle.  The  grid 
in  this  procedure  is  larger  than  in  the  conventional  one  (l6  x  48  instead  of 
10  x  30),  requiring  a  somewhat  larger  core,  but  the  computing  time  for  finding 
a  slip  line  field  that  meets  the  boundary  condition  is  much  less. 

The  interface  friction  angle,  as  shown  in  the  flow  diagram,  is  assumed 
to  be  constant.  However,  the  program  can  as  well  accommodate  a  6  angle  that 
varies  along  the  interface. 

C.  Computation  of  a  Matching  Set  of  Slip  Line  Fields 

For  the  problem  of  soil-wheel  interaction,  a  matching  set  of  slip  line 
fields  must  be  found  that  meets  the  constraints  described  in  Section  II. C. 

For  this  purpose,  it  is  convenient  to  start  the  computations  with  the  rear 
slip  line  field  using  the  subroutine  outlined  in  Section  III.B,  but  allowing 
for  appropriate  sign  changes  due  to  the  fact  that  the  rear  field  is  a  mirror 
image  of  the  front  field.  The  reason  for  starting  the  computation  with  the 
rear  field  is  twofold:  first,  the  rear  angle  varies  within  a  narrow  range; 
second,  a  normal  stress  from  the  rear  field  (q^)  at  O’  can  be  generally 
matched  by  a  normal  srress  from  the  forward  field.  '  This  is  not  true  vice- 
versa.  When  q^  for  the  assumed  rear  angle  (or  )  and  interface  friction  angle 
(6)  is  determined,  the  front  field  is  found  by  varying  the  entry  angle  (or) 
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until  the  normal  stress  from  the  front  field  (q^)  matches  q^  within  the 
allowed  tolerance.  Since  the  interface  friction  angle,  6,  is  assumed  to  be 
the  same  for  both  the  forward  and  rear  fields,  the  shear  stress  at  o/  from 
both  fields  is  the  same  (within  the  tolerance)  when  q^  and  q^  are  matched. 
The  flow  diagram  for  this  procedure  is  shown  in  Fig.  When  a  matching  set 
of  rear  and  forward  fields  is  found  by  this  procedure,  the  load,  torque, 
and  drawbar  pull  are  computed  by  appropriate  numerical  integration  of  the 
interface  stresses  for  the  assumed  value  of  a  and  6.  The  slip  is  computed 
from  Eq.  (13)  and  the  sihkage  from  the  entry  angle. 


D.  Inversion  Procedure 

The  procedure  described  in  Section  III.C  yields  the  wheel  performance 
parameters  for  an  assumed  rear  angle  and  interface  friction.  An  inversion 
procedure  is  required  to  solve  for  torque,  slip,  and  sinkage  when  wheel  load, 
drawbar  pull,  and  soil  properties  are  given.  Such  a  procedure  is  outlined 
below.  Load,  drawbar  pull,  and  torque  are  functions  of  the  rear  angle,  or  , 
and  interface  friction  angle,  6,  as  expressed  by  the  following  relationships: 

L  =  f ±  (Q'r,  6) 


DB  -  f2  (orr,  6)  (19) 

T  =  f 3  (*.,  6) 

The  functions  f^,  fg,  f  are,  of  course,  not  closed  form  functions.  In  the 

solution  procedure  outlined  in  Section  III. B,  one  set  of  L,  DB  and  T  values 

are  found  for  an  assumed  or  and  slip  corresponding  to  an  assumed  6.  The 

procedure  to  find  the  torque  for  given  load  and  drawbar  pull  consists  of 

finding  a  and  6  which  yield  the  given  load  and  drawbar  pull;  once  the 
r  > 

matching  set  of  slip  line  fields  for  these  conditions  is  found,  the  torque, 
slip,  and  sinkage  is  also  available  from  the  computations.  Even  if  the 
functions  f^  and  f  were  known,  the  solution  for  L  and  DB  would  require  the 
solution  of  a  system  of  two  nonlinear  equations.  Since  there  is  no  generally 
valid  theorem  for  the  solution  of  this  problem,  and  convergence  criteria  for 
iterative  solutions  cannot  be  established  if  the  derivatives  of  f^  and  fg  are 
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not  known,  it  was  necessary  to  study  the  general  behavior  of  these  functions 
and  to  make  judicious  use  of  some  of  their  properties  to  devise  an  efficient 
and  convergent  iteration  scheme  for  the  solution  of  the  problem. 

There  are  two  properties  of  the  wheel  performance  relationships  ex¬ 
pressed  by  Eqs.  (19)  that  are  useful  for  the  solution  of  the  inversion 
problem.  First,  the  wheel  load  increases  monotonically  with  the  rear  angle 
for  a  constant  angle  of  interface  friction.  Second,  the  pull  coefficient 
DB/L,  often  used  as  a  dimensionless  parameter  in  wheel  performance  studies, 
was  found  to  increase  monotonically  with  the  interface  friction  angle.  Al¬ 
though  there  is  no  rigorous  proof  for  this  to  hold  true  for  every  combination 
of  conditions,  theoretical  considerations  and  evaluation  of  computations 
performed  for  a  wide  range  of  conditions  indicate  that  this  second  relation¬ 
ship  is  generally  valid.  In  the  inversion  procedure,  these  two  relationships 
constitute  the  basis  of  the  iteration  procedure,  which  is  shown  schematically 
in  Fig.  10.  The  iteration  steps  are  described  in  greater  detail  in  Section 
III.F  where  the  comprehensive  computer  program  is  described. 

E.  Computation  of  Wheel  Performance  in  Purely  Cohesive  Soils 

The  computational  procedures  described  in  the  preceding  paragraphs  refer 
to  wheel  performance  calculations  where  the  numerical  solution  of  differential 
Eqs.  (3)  is  required.  As  it  was  pointed  out  in  Section  II. D,  the  form  of  the 
governing  differential  equations  (Eqs.  4)  is  somewhat  different  in  purely  co¬ 
hesive  soils.  The  solution  of  these  equations  yields  the  slip  line  field 
geometry  and  associated  stresses,  just  as  Eqs.  (3)  do  in  the  case  of  frictional 
soils.  A  study  of  the  slip  line  fields  and  interface  stresses  obtained  from 
the  numerical  solution  of  Eqs.  (4)  showed,  however,  that  it  is  not  necessary 
in  this  special  case  to  perform  these  numerical  calculations  because  the  inter¬ 
face  stresses  can  be  calculated  with  very  close  approximation  from  the  following 
formulas  (Ref.  3 ) : 

a  =  c  (l  +  tt  +  cos  6  ±  ^  6  +  o') 

n  c  2  c 

(20) 

t  =  c  sin  6 

c 

(The  upper  sign  refers  to  the  front,  and  the  lower  to  the  rear  field). 
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In  the  above  equations,  the  angle  6^  is  equivalent  to  the  interface 
friction  angle.  It  is  defined  as 

c  .  Tmob  \ 

6  =  arc  sin  -  (21) 

c  T 

max 

The  requirement  that  the  normal  stress  from  the  forward  and  the  rear  field 
be  equal  at  the  angle  of  separation  leads  to  the  following  relationship: 

cv  =  6  ,  (22) 

m  c/2 

This  relationship  is  not  consistent  with  Eq.  (l6),  established  for  c  -  cp 
soils,  and  this  apparent  contradiction  is  the  consequence  of  assuming  uniform 
interface  friction  along  the  interface.  In  actuality,  the  interface  friction 
is  not  uniform  and  adjusts  itself  in  such  a  way  that  the  two  conditions  are 
met  at  the  angle  of  separation.  Therefore,  the  assumption  of  a  uniformly 
developed  interface  friction  is  retained  for  computational  purposes. 

In  contrast  to  frictional  soils,  in  the  case  of  cp  =  0,  the  interface 
stresses  at  any  point  depend  only  on  the  central  angle  and  interface  friction 
but  not  on  the  entry  or  rear  angle.  For  this  reason,  further  assumption  re¬ 
garding  these  angles  is  necessary  to  make  the  problem  definite.  It  is  assumed 
on  the  basis  of  experiments  that  the  rear  angle  is  one  third  of  the  entry 
angle.  The  computation  of  wheel  performance  parameters,  based  on  this  assump¬ 
tion,  follows  the  flow  chart  shown  in  Fig.  11. 

F.  Description  of  Comprehensive  Computer  Program 

In  the  following  description,  emphasis  is  laid  on  the  user's  anticipated 
needs  rather  than  on  details  of  the  program.  The  comprehensive  program  con¬ 
sists  of  the  following  programs: 

Main  program  LKWH 

Subroutine  1  SLIP 

Subroutine  2  LKWC 

Listings  of  the  above  programs  are  given  in  the  Appendix. 
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1.  Subscripted  Variables  and  Dimension  Statements 
Main  program  LKWH: 

J  designates  a  location  at  the  interface. 

K  designates  the  numeral  assigned  to  consecutive  iterations  in  the 
inversion  procedure  Section  III.D. 

HH(J)  =  a 

QQ(J)  =  q, 

J 

EE(J)  =  T. 

3 

WE(j),  DRB(j),  TRQ(j)  =  auxiliary  variables 

LD(K)  =  Load 

DR(K)  =  Drawbar  pull 

TR(K)  =  Torque 

AR(k)  =  Rear  angle 

AE(K)  =  Entry  angle 

DEL(k)  =  Interface  friction  angle 

PU(k)  =  Pull  coefficient 

LC(K)  =  Load  coefficient 

Subroutine  SLIP: 

The  variables  in  Eqs.  (l6)  are  designated  as  follows: 

X(I,J)  =  x.  . 

Z(I,J)  =  z.  . 

1.3 

S(I,J)  =  o.  . 

1.3 

T(I,J)  =  6  , 

The  dimension  statement  corresponds  to  a  48  x  16  grid  for  the  computation 
of  the  geometry  of  slip  line  fields.  All  J  locations  are  not  necessarily  used 
in  actual  computations,  as  indicated  in  Section  III.B.  Several  hundred  com¬ 
putations  were  performed  satisfactorily  with  the  above  grid  size  and,  therefore, 
no  need  for  the  change  of  the  grid  size  is  anticipated. 
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Variables  with  one  subscript  are  as  follows: 

D(j)  -  6. 

a 

h(j)  =  o', 

o 

Q(J)  =  q. 

J 

E(j )  =  T. 

3 

Where  j  designates  a  location  at  the  interface.  A(j),  B(j),  C(j)  are  auxiliary- 
variables. 

The  dimensions  of  the  variables  with  one  subscript  are  tied  to  the  J 
dimensions  of  the  variables  with  two  subscripts. 

Subroutine  LKWC: 

Designation  of  subscripted  variables  are  the  same  as  in  the  main  program. 

2.  Input  Files 

Input  variables  and  constants  that  the  user  may  want  to  change  are  read 
from  input  files  so  that  the  program  need  not  be  recompiled  if  these  data  are 
changed.  Three  input  files  are  used  in  the  present  program;  they  may  be  com¬ 
bined  in  one  if  so  desired.  The  input  files  contain  the  following  data: 

Input  file  SOL  contains  data  on  soil  properties  in  the  following  order: 

Cohesion  (CO)  in  lbs/sq  ft 
Friction  angle  (FO)  in  degrees 
Unit  weight  of  soil  (GO)  in  lbs/cu  ft 
Slope  angle  (SO)  in  degrees 

Slip  parameter  j  (SU),  dimensionless  number 
Slip  parameter  K  (SLK),  dimensionless  number 

It  is  noted  here  that  cohesion  and  friction  angle  are  strength  parameters 
as  determined  by  triaxial  tests  at  a  rate  of  loading  comparable  to  that  ob¬ 
tained  under  a  moving  wheel.  Strength  parameters  from  direct  shear  tests  may 
also  be  used.  Strength  parameters  obtained  from  Bevameter  shear  tests, 
sheargraph,  and  other  devices  do  not  yield  true  internal  strength  parameters 
and  are  not  to  be  used  with  this  program. 
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Input  file  WH1  contains  data  on  wheel  geometry  and  wheel  forces  in  the 
following  order: 

Wheel  radius  (RO)  in  ft 
Wheel  width  (BO)  in  ft 
Wheel  load  (LO)  in  lbs 
Drawbar  pull  (DB)  in  lbs 

Input  file  TOL  contains  tolerance  limits  and  numerical  constants  that 
determine  the  magnitude  of  changes  in  the  rear  and  interface  friction  angle 
in  the  inversion  procedure  and  the  maximum  number  of  iterations  with  these 
constants.  The  values  of  limits  and  constants  in  this  input  file  have  been 
selected  on  the  basis  of  experience  with  the  program  and  they  may  be  changed 
judiciously  according  to  required  accuracy,  allotted  computer  time,  etc. 

Following  is  a  list  of  data  in  this  file  in  the  order  they  are  read  with 


recommended  values : 

Tolerance  limit  for  load  (TOL)  0.1 

Tolerance  limit  for  pull  coefficient  (TOP)  0.03 

Numerical  constants  (Bl)  1.09 

for  the  determination  (B2)  0.95 

of  the  size  of  6  and  a  (b4)  20.0 

increment  s  ( B5 )  2.0 

(B6)  2.0 

(B7)  0.9 

(B8)  0.7 


Limit  for  the  number  of  iterations  1st  Phase  (K3)  16 

2nd  Phase  (K5)  12 

The  tolerance  limits  define  the  range  of  accuracy  desired  in  the  inversion 
procedure.  There  are  inherent  inaccuracies  in  prediction  theories  due  to  the 
various  assumptions,  inaccuracies  in  the  determination  of  soil  properties  and 
the  approximations  inherent  in  numerical  methods.  Thus,  there  is  no  point  in 
performing  iterations  to  achieve  an  apparent  accuracy  in  the  prediction  when 
inaccuracies  from  other  sources  would  dominate. 
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The  recommended  tolerance  limit  for  load  (TOL)  is  ±  10  percent,  and  for 
the  pull  coefficient  (TOP)  ±  3  percent*  The  numerical  constants  listed  are 
consistent  with  these  limits.  Their  significance  is  explained  in  the  descrip¬ 
tion  of  the  inversion  procedure.  Because  of  the  nature  of  the  problem  of  the 
inversion  procedure,  it  was  not  possible  to  establish  criteria  for  the  con¬ 
vergence  of  the  iteration  procedure.  Instead,  a  limit  (K3)  is  set  for  the 
number  of  iterations  performed  in  steps  as  determined  by  one  set  of  constants 
set  forth  in  the  input  file.  If  solution  is  not  reached  within  the  allotted 
number  of  iterations,  the  program  changes  these  constants  and  performs  a 
maximum  of  K5  additional  iterations  where  the  steps  correspond  to  the  changed 
constants.  The  recommended  limits  for  iterations  may  be  changed  but  only 
within  the  limits  of  dimension  statements. 

3.  Decisions,  Assumptions  and  Iterations  in  the  Main  Program 

The  general  logic  of  the  main  program  is  shown  in  the  flow  diagrams . 

There  are,  however,  some  provisions  in  the  main  program  that,  for  clarity, 
have  not  been  included  in  the  flow  diagrams.  These  are  discussed  below. 

In  Section  III.E,  the  special  case  of  cp  =  0  was  discussed.  It  was  found 
that  the  simplified  method  applicable  to  this  case  also  yielded  reasonably 
accurate  results  for  small  cp  angles  if  an  appropriate  correction  in  the  value 
of  cohesion  was  made.  Therefore,  a  provision  was  made  in  the  main  program  to 
perform  calculations  by  the  subroutine  LKWC  whenever  cp  <12  degrees. 

It  is  also  assumed  that  soils  that  exhibit  a  friction  angle  greater  than 
12  degrees  are  only  partially  saturated  and,  therefore,  are  liable  to  change 
their  unit  volume  (void  ratio  or  porosity)  under  the  action  of  the  wheel. 

While  12  degrees  were  selected  on  the  basis  of  judgment,  it  is  obviously  an 
arbitrary  average  value  introduced  in  the  computer  program  as  a  necessity.  The 
consequence  of  the  volume  change  of  soil  under  the  action  of  the  wheel  is  that 
the  soil  in  the  rear  field  is  precompressed  by  the  stresses  in  the  front  field; 
therefore,  generally  the  strength  of  soil  is  higher  in  the  rear  field  than  in 
the  front. 

It  is  desirable  for  prediction  purposes  to  limit  soil  characterization 
to  the  strength  parameters  of  the  undisturbed  ground.  To  improve  the  accuracy 
of  predictions,  however,  it  is  necessary  to  consider  the  change  in  the  strength 
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properties  due  to  the  compacting  effect  of  the  wheel.  To  reconcile  these 
conflicting  requirements,  a  crude  estimate  of  the  change  in  strength  proper¬ 
ties  is  made  according  to  the  following  scheme: 


Strength 

Parameters  of 

Qualitative 

Undistrubed 

Assessment 

Ground 

of  Soil 

(A) 

0  >  ho 

dense 

C  <  50  lbs/sq.ft. 

frictional 

(B) 

0  <  h0 

loose  and 

C  <  50  lbs/sq.ft. 

medium  dense 
frictional 

(c) 

C  >  50  lbs/sq.ft. 

cohesive  soil 

Effect  of  Wheel  Passage  on 
Soil  Strength  in  Rear  field 


no  change 

friction  angle  and  High  sinkage: 
unit  weight  increase  increase  10$ 
due  to  densification  Low  sinkage: 

increase  5$ 

effect  of  wheel  passage  similar  to 
that  of  overconsolidation, 
cohesion  increases,  friction  angle 
decreases 


While  the  determination  of  whether  the  soil  is  in  the  A,  B,  or  C  category 
is  straightforward,  further  considerations  are  necessary  to  estimate  the  sink- 
age  beforehand  in  category  B  and  the  strength  properties  in  category  C. 

To  this  end,  a  measure  of  the  maximum  vertical  pressure  (SIM)  is  estab¬ 
lished  by  dividing  the  wheel  load  with  a  contact  area  corresponding  to  0.8 
arc  length.  This  is  compared  to  an  estimate  of  the  maximum  normal  stress 
obtained  from  the  following  formula: 

CTn  =  (C+0.33  VR)  Nq6  (23) 


where  N 


q6 


is  a  bearing  capacity  factor  computed  from  the  following  fomula: 


N  R  = 
q6 


cos  6  +  cos  6 


2 

cos  cp 


2  2 

cos  6  -  cos  6  +  cos  cp 


2  tancp  (9p-V  (24) 


If  the  estimated  maximum  normal  stress  is  lower  than  the  estimated  maxi¬ 
mum  vertical  pressure,  low  sinkage  is  anticipated  and  a  5  percent  increase  is 
indicated.  Otherwise,  the  increase  is  10  percent.  If  the  strength  parameters 
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of  the  undisturbed  ground  fall  in  category  C,  then  the  strength  parameters 
in  the  rear  field  are  assumed  as  shown  in  Fig.  12.  In  this  case,  the  strength 
of  soil  in  the  rear  field  is  higher  than  in  the  front  for  normal  pressures 
lower  than  the  estimated  maximum  normal  stress. 

The  inversion  procedure  discussed  in  Section  III.  D  requires  an  estimate 
of  initial  data  for  the  first  run.  These  were  chosen  as  ot ^  =  10  degrees  and 

6  =  (0.35  +  (DB/L)).  cp  (25) 

In  the  inversion  procedure,  the  values  of  and  6  are  changed  as  shown 
in  the  flow  diagram  in  Fig.  10.  The  magnitude  of  the  changes  in  each  step 
is  determined  by  a  formula  of  the  general  type,  as  follows: 

Vi E  CA  -  c2  -  V  +  c3 

(26) 

“  Cl"x  +  W  +  C5 

The  constants  in  the  above  equations  are  chosen  appropriately  to  result  in 
a  change  indicated  in  the  flow  chart  in  Fig.  10.  The  following  list  gives 
the  options  for  the  above  constants: 

C1  =  B1  ,  B2 

C2  =  B3 

c3  =  B4  ,  B8 

%  =  B5 

C5  =  b6  ,  B8 


The  following  absolute  limits  have  been  set  in  the  computer  program  for 


the  variables  ot  , 
r’ 


and  5  : 
e 

O' 

r  max 

(AMAX)  = 

c 

H 

CM 

O' 

r  mm 

(AMIN)  = 

i° 

ot 

e  max 

(AEMAX)  = 

=  6o° 

6 

max 

(DMAX)  = 

O.98  x  arc  tan  (sin  0) 

6  . 

mm 

(EMIN)  = 

0.05  0 
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The  limits  on  the  entry  and  rear  angles  are  based  on  available  experi¬ 
mental  data,  and  are  self-explanatory.  A  factor  of  O.98  is  applied  to  the 
theoretical  6  value  to  avoid  certain  computations  close  to  limit  conditions 

where  formulas  become  indefinite.  The  6  .  limit  is  arbitrary  and  serves  to 
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avoid  computations  that  would  not  be  useful  in  approaching  the  solution. 

The  above  limits  are  also  used  as  criteria  for  an  acceptable  solution. 

If,  for  any  particular  input  condition,  the  flow  chart  in  Fig.  10  indicates 
a  need  for  a  change  in  either  of  the  variables  6  and  a  ,  but  it  is  not  possible 
to  execute  that  change  because  of  these  limits,  then  there  is  no  solution  for 
the  input  condition.  If  the  iteration  procedure  is  terminated  because  the 
maximum  limits  would  have  to  be  exceeded  for  a  solution,  the  "no  solution" 
condition  indicates  a  "no  go"  condition,  or  100  percent  slip.  If  the  iteration 
procedure  is  terminated  because  the  minimum  limits  would  have  to  be  transgressed, 
the  "no  solution"  condition  indicates  that  the  soil  strength  is  so  high  that 
failure  in  the  soil  does  not  develop  under  the  input  load  and  a  "hard  soil" 
condition  exists. 

The  results  of  each  step  in  the  computer  program  iteration  are  preserved 
as  subscripted  variables  with  the  subscripts  denoting  the  numerals  assigned 
to  the  consecutive  iterations.  If  the  solution  is  not  reached  within  K3 
iterations,  the  program  changes  the  steps  so  that  the  increments  or  decrements 
in  the  6  and  01^  values  become  about  half  of  that  during  the  first  K3  iterations 
and  performs  an  additional  K5  iterations  with  changed  constants  that  reduce  the 
magnitude  of  6  and  or  steps  to  about  50  percent  of  that  in  the  first  K3  itera¬ 
tions.  If  no  solution  is  reached  within  K3  and  K5  iterations,  the  program 
scans  the  results  and  computes  the  output  data  by  an  approximate  interpolation, 
if  appropriate. 

Output  -  Output  data  in  the  program  are  in  printout  form  that  can  be 
easily  changed  to  other  output  formats  if  desired.  Following  are  the  output 
data  and  symbols: 


Load 

LD(K) 

Drawbar  pull 

DR(K) 

Torque 

TR(K) 

Slip 

SLP 

Sinkage 

SNK 

23 


The  load  and  drawbar  pull  values  are  printed  out  for  information  only 
since  they  may  differ  from  the  input  load  and  drawbar  pull  values  by  the 
tolerance  limits.  The  torque  and  slip  values  may  be  used  in  conjunction 
with  the  transmission  and  engine  subroutines  available  at  TACOM  (Ref.  9)  to 
determine  vehicle  speed,  power  requirements,  and  fuel  consumption. 
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IV.  EXPERIMENTAL  PROGRAM 


A.  Performance  of  Experiments  at  the  Davidson  Laboratory  of  Stevens  Institute 
of  Technology 

1.  Test  Facility  and  Equipment 

a.  Soil  Bins 

The  soil  bin  of  the  Davidson  Laboratory  is  4o  feet  long  and  7  feet  wide. 

The  bin  is  divided  in  half  longitudinally  to  form  two  separate  40'  x  3  1/2'  bins. 
One  bin  was  filled  with  sand,  the  other  with  loam,  both  to  a  depth  of  approxi¬ 
mately  21/2  feet. 

b.  Dynamometer  Carriage 

The  dynamometer  carriage  contains  the  test  wheel,  the  wheel  drive  motor, 
the  hydraulics  required  to  drive  the  wheel,  the  wheel  loading  system,  and  the 
dynamometer  balance.  The  dynamometer  carriage  is  unpowered;  carriage  motion  is 
supplied  by  connecting  it  to  an  auxiliary  carriage  that  contains  the  soil 
processing  equipment  and  that  is  propelled  by  an  off-carriage  chain  drive  system. 

c.  Test  Wheel 

The  wheel  used  was  made  of  plywood,  28  inches  in  diameter  and  4  1/2  inches 
wide.  The  wheel  was  attached  to  the  wheel  drive  motor  by  a  metal  faceplate.  On 
the  circumference  of  the  wheel  were  mounted  four  sensors,  each  to  measure  the 
normal  and  tangential  forces  at  the  periphery  of  the  wheel.  These  sensors  were 
cantilevered,  L-shaped  arms  (Figure  13).  Each  leg  of  the  L  had  been  strain- 
gauged  and  wired  in  such  a  way  that  the  leg  perpendicular  to  the  circumference 
measured  the  tangential  forces;  the  other  leg  measured  the  normal  forces.  The 
four  sensors  were  mounted  transversely  across  the  face  of  the  wheel  from  the 
center  to  one  edge  to  give  a  record  of  the  cross-wise  distribution  of  the  load 
on  the  wheel  as  it  rolls  through  the  soil.  The  signals  from  the  sensors  were 
transmitted  through  a  slip  ring  mounted  on  the  axle  to  overhead  cables  and  then 
to  the  recorders. 

The  head  of  each  sensor  was  a  square  plate  approximately  3/8-inch  on  a  side. 
Each  head  was  tangent  to  the  circumference  of  the  wheel  and  protruded  through 
the  approximate  center  of  a  l/2-inch  metal  grid.  Details  of  the  sensor  con¬ 
figuration  within  the  grid  and  the  spacing  of  the  sensors  on  the  wheel  face  are 
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shown  in  Figure  l4.  To  prevent  soil  from  getting  between  the  sensors  and  the 
grid,  a  thin  rubber  membrane  was  cemented  around  the  circumference. 

d.  Wheel  Drive 

The  wheel  was  driven  directly  from  a  hydraulic  motor  that  was,  in  turn, 
connected  to  an  electrically-driven,  variable -displacement  hydraulic  pump. 
Varying  the  pump  output,  therefore,  varied  the  wheel  speed  in  direct  proportion. 
Difficulty  was  encountered  in  maintaining  a  constantly  smooth  rotation  of  the 
wheel  due  to  a  porting  action  at  the  motor.  This  is  attributed  to  the  fact  that 
the  wheel  pump  was  designed  to  operate  at  speeds  considerably  greater  than  those 
used  in  the  testing  program. 

e.  Dynamometer 

The  wheel  motor  was  mounted  directly  onto  a  six-component  dynamometer  that 
measured  all  three  forces  (load,  drawbar  pull,  and  side-force)  all  three  moments 
(yaw  moment,  input  torque,  and  roll  moment). 

f.  Wheel  Loading  System 

The  dynamometer  was  mounted  directly  to  a  loading  device.  Loads  on  the 
wheel  were  applied  by  an  air-actuated  belofram  and  servo  system.  Air  pressure 
above  and  below  the  belofram  controlled  the  load.  As  the  wheel  sank  into  the 
soil,  a  hydraulic  servo  system  positioned  the  beloframs  in  response  to  this 
wheel  movement  so  that  the  beloframs  were  always  near  to  their  center  position, 
regardless  of  the  sinkage  experienced.  The  entire  wheel  and  loading  system 
could  be  raised  or  lowered  by  manually  actuating  the  servo  system,  thus  clearing 
the  wheel  from  the  soil  for  soil  processing  or  wheel  maintenance. 

g.  Wheel  Speed  and  Sensor  Position 

Wheel  speed  was  measured  in  two  ways.  A  tachometer  connected  directly  to 
the  wheel  drive  motor  indicated  the  wheel  speed.  Wheel  speed  could  also  be 
calculated  from  angular  position  indicators  mounted  on  the  wheel.  These  position 
indicators,  mounted  every  two  degrees  about  the  wheel  generated  spikes  on  the 
recording  traces.  These  spikes  began  when  the  sensors  were  50°  before  bottom 
center  and  ended  when  they  were  50°  after  bottom  center. 

Initially,  angular  position  was  indicated  by  a  series  of  brass  screws 
placed  at  2°  intervals  in  a  100°  arc.  Two  separate  bronze  spring  brushes  were 
placed  side  by  side  so  that  both  brushes  touched  each  screw  as  it  passed.  This 
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closed  a  circuit  that  generated  a  spike  on  both  strip  chart  records. 

The  speed  of  the  wheel  frequently  caused  the  springs  to  bounce,  and  to 
miss  a  few  spikes.  Therefore,  after  the  sand  tests,  a  new  position-indicating 
device  was  employed.  A  semicircular  piece  of  sheet  metal  was  slit  every  2°  and 
was  mounted  to  the  side  of  the  wheel.  A  photocell  detected  these  slits  as  they 
passed  by,  and  caused  a  similar  spike  to  be  recorded. 

The  spikes  thus  generated  could  then  be  used  to  locate  accurately  the 
position  of  the  sensors  relative  to  the  soil  and,  by  measuring  the  distance 
between  the  spikes  and  the  chart  paper  speed,  to  compute  the  wheel  speed  while 
the  sensors  were  passing  through  the  soil. 

h.  Carriage  Velocity 

The  carriage  velocity  was  measured  by  signals  from  markers  spaced  one 
foot  apart  along  the  side  of  the  soil  bin.  As  the  carriage  travels  down  the  bin, 
the  markers  trip  a  microswitch  that,  in  turn,  causes  a  spike  to  appear  on  the 
strip  chart  record.  With  this  record  appearing  every  foot  of  carriage  travel, 
and  a  known  paper  speed,  the  carriage  velocity  could  be  easily  calculated.  For 
all  tests,  carriage  speed  was  held  constant;  slip  was  controlled  by  changing 
wheel  speed. 

i .  Sinkage 

The  sinkage,  or  vertical  travel  relative  to  the  soil  surface,  was  measured 
by  a  multiple  turn  potentiometer.  A  spring  was  secured  at  one  end  of  the 
carriage  and  then  wrapped  around  a  pulley  on  the  potentiometer,  which  was 
fastened  to  the  wheel  mounting  apparatus.  When  the  wheel  moved  in  the  vertical 
direction,  the  string  caused  the  pulley  to  rotate,  thereby  giving  a  signal  that 
was  calibrated  to  wheel  sinkage.  Zero  sinkage  is  established  with  the  bottom  of 
the  wheel,  just  touching  the  undisturbed  soil  surface. 

j .  Instrumentation 

The  instruments  used  were  Sanborn  Models  150  and  85O  multi-channel  strip 
chart  recorders.  The  8-sensor  signals  were  fed  into  the  Sanborn  85O;  the  Sanborn 
150  recorded  wheel  speed,  carriage  speed,  sinkage,  torque,  horizontal  drawbar 
pull,  and  vertical  load.  Both  recorders  registered  the  sensor  position  spikes. 
Before  each  day's  testing,  a  calibration  was  made  of  all  sensors  so  that  the 
strip  charts  produced  a  signal  accurate  to  one  part  in  forty. 
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2.  Calibration  Technique 

a.  Calibration  of  Sensors 

Calibration  of  the  normal  load  on  the  sensors  was  done  by  balancing  1-  to 
5-pound  loads  on  the  face  of  the  sensor.  A  small  nut  was  placed  on  top  of  the 
sensor  while  the  sensor  was  in  the  top  center  position,  and  then  weights  were 
balanced  on  the  nut  to  prevent  contact  of  the  weight  with  the  wheel  surface. 

This  process  was  done  for  each  of  the  four  sensors. 

Calibration  of  the  tangential  loads  was  performed  by  positioning  the  wheel 
so  that  the  sensors  were  at  90°  from  the  vertical.  While  in  this  position,  a 
small  pin  was  placed  in  a  hole  in  each  sensor.  Then  a  string  carrying  a  1-, 

2-,  or  3-pound  weight  was  hung  on  that  pin. 

b.  Calibration  of  Dynamometer 

Calibration  of  the  drawbar  pull  force  was  done  by  placing  a  harness  on 
the  'wheel  and  leading  a  wire  cable  horizontally  over  a  pulley.  Weights  were 
then  hung  on  this  cable.  Vertical  load  was  obtained  by  placing  the  wheel  on  a 
pre -calibrated  load  cell  and  loading  the  beloframs  with  air  with  the  servos  in 
the  automatic  mode.  Since  it  was  difficult  to  maintain  a  steady  load  with  the 
equipment,  load  calibration  was  made  before  and  after  every  test.  Torque  was 
calibrated  by  setting  the  load  cell  one  foot  from  the  axle  center  under  one  arm 
of  a  beam  that  was  bolted  to  the  axle  of  the  wheel.  The  sihkage  was  calibrated 
by  measuring  the  wheel  displacement  with  a  steel  rule. 

3.  Preparation  of  Soil  Bed 

The  following  three  types  of  soil  beds  were  prepared: 

Loose  sand 
Danse  sand 
Loam 

The  loose  sand  was  processed  by  a  gyrotiller,  a  leveling  blade,  and,  for 
dense  sand  tests,  a  plate  vibrator.  The  gyrotiller  disturbed  the  sand  to  a 
depth  of  18  inches  and  produced  an  uniformly  loose  sand  layer.  The  plate 
vibrator  was  used  for  compacting  the  sand;  it  was  towed  at  a  constant  speed 
behind  the  carriage  after  tilling.  It  compacted  the  sand  to  about  a  6-inch  depth. 


28 


The  loam  bed  was  prepared  at  a  moisture  content  of  approximately  16  percent. 
The  loam  was  processed  by  the  gyrotiller,  leveling  blade,  and  a  roller.  First, 
the  gyrotiller  tilled  the  clay  to  a  depth  of  six  inches  while  the  blade  leveled 
and  smooth  the  soil.  Next,  a  lawn  roller  filled  with  water  to  about  two-thirds 
capacity  was  used  to  compact  the  loam  (Figure  15).  The  loam  bed  was  protected 
from  evaporation  by  a  plastic  cover;  water  was  added  by  sprinkling  whenever  loss 
of  moisture  was  observed. 

k.  Control  Tests  in  Soil  Bed 

Cone  penetrometer  tests  were  performed  in  the  soil  bed  to  check  the  uni¬ 
formity  and  condition  of  the  soil  before  the  wheel  performance  tests  (Figure  16). 
Cone  penetrometer  tests  were  also  performed  for  selected  tests  in  the  ruts  after 
the  passage  of  the  wheel.  Results  of  these  tests  are  reported  in  Section  V. 

The  moisture  content  of  the  load  bed  was  determined  during  preparation  at 
several  locations  and  undisturbed  cylinder  samples  were  taken  to  determine  the 
density  of  the  soil  bed.  These  data  were  essential  for  the  proper  duplication 
of  soil  conditions  in  the  triaxial  tests  performed  for  the  determination  of  shear 
strength  properties. 

5.  Strength  Properties  of  Test  Soils 

For  the  determination  of  strength  properties  of  the  test  soils,  triaxial 
tests  were  performed  on  them  in  the  Soil  Mechanics  Laboratory  at  Grumman. 

Results  of  these  tests  are  summarized  below. 

a.  Sand 

The  strength  properties  of  sand  were  determined  by  two  series  of  triaxial 
tests,  one  performed  at  a  low  density  (98  lbs/cu.  ft.)  and  the  other  one  at  a 
high  density  (106  lbs/cu.  ft.).  The  air  dry  sand  was  placed  in  one-inch  layers 
in  a  7-inch-high,  2.8-inch-diameter  mold  to  ensure  uniform  density  throughout 
the  sample.  During  the  tests,  the  samples  were  allowed  to  change  their  volume 
freely;  volume  changes  were  computed  from  circumferential  gauge  readings.  The 
tests  were  stress  controlled;  load  increments  were  applied  after  the  stabilization 
of  vertical  displacement.  The  sand  was  found  to  be  insensitive  to  the  rate  of 
loading;  therefore,  no  attempt  was  made  to  duplicate  the  loading  rates  in  the 
wheel  tests.  The  sand  failed  at  relatively  low  vertical  strain  (less  than  5  per¬ 
cent  in  each  test).  Mohr  circle  representation  of  the  triaxial  test  results  is 
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shown  in  Figures  17  and  18.  It  is  seen  that  the  strength  envelope  is  slightly- 
curved  and  the  friction  angle  at  low  normal  stresses  is  relatively  high. 

Grain  size  distribution  of  the  sand  is  shown  in  Figure  19. 

b.  Loam 

The  moisture  samples  taken  from  the  test  bed  showed  a  variation  of 
moisture  content  from  14.7  to  16.4  percent.  The  dry  density  of  the  loam, 
determined  from  cylinder  samples,  varied  from  I.29  to  1.33  g/cu.  cm.  For  the 
determination  of  the  strength  properties  of  the  soil  bed,  triaxial  test  samples 
were  prepared  in  a  7-inch-high,  2.8-inch-diameter  mold  in  1-inch-thick  layers 
from  uniformly  mixed  soil  kept  at  16  percent  moisture  content.  Some  moisture 
content  was  lost  during  the  preparation  and  the  actual  moisture  content  of  the 
samples  was  somewhat  lower.  Preliminary  tests  indicated  that  the  loam  at  this 
moisture  content  was  moderately  sensitive  to  the  rate  of  loading;  the  final  tests, 
therefore,  were  performed  as  rapidly  as  it  was  possible  with  the  available  stress 
controlled  triaxial  apparatus.  Failure  of  the  loam,  as  expected,  occurred  at 
relatively  high  strains,  as  shown  in  Figure  20.  Mohr  circle  representation  of 
the  triaxial  test  results  is  shown  in  Figure  21.  It  is  noted  that  the  moisture 
content  of  Test  Number  16  was  slightly  higher  and  its  dry  density  slightly 
lower  than  that  of  the  other  tests  in  the  series,  resulting  in  a  relatively  low 
strength  that  was  not  considered  in  drawing  the  strength  envelope. 

c.  Laboratory  Cone  Penetrometer  Tests 

To  correlate  the  strength  tests  with  the  cone  penetrometer  tests  performed 
in  the  soil  bed,  cone  penetrometer  tests  were  performed  in  laboratory  soil  bins 
under  controlled  conditions.  Results  of  these  tests  for  sand  are  shown  in 
Figure  22.  These  tests  show  that  penetration  resistance  increases  with  depth  as 
expected.  While  the  increase  is  approximately  linear  at  the  lower  densities  it 
is  not  possible  at  higher  densities  to  use  linear  approximation  for  the  cone 
penetration  curve.  Nevertheless,  these  graphs  can  be  used  as  guides  to  estimate 
the  in  situ  density  of  soil  bed  and,  therefrom,  its  strength  properties. 

A  typical  results  of  a  cone  penetrometer  test  in  the  loam  is  shown  in 
Figure  23.  The  cone  penetration  resistance  reaches  an  approximately  constant 
value  after  about  two  inches  of  penetration.  Cone  penetrometer  tests  were  per¬ 
formed  in  a  laboratory  bin  where  the  loam  was  prepared  at  about  15  percent 
moisture  content  at  various  densities.  The  constant  value  of  cone  penetration 
resistance  reached  after  an  initial  increasing  portion  is  plotted  in  Figure  24. 
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6.  Performance  of  Test 

a.  Test  Plan 

The  test  plan  was  to  obtain  three  different  slip  ranges  for  three  different 
load  ranges  for  each  of  three  different  soil  conditions.  The  slip  and  load 
ranges  were  as  follows:  > 

Slip  Ranges  for  all  Tests 
Low  5 $  -  10$ 

Medium  15$  -  25 $ 

High  greater  than  25$ 

Load  Ranges  for  Loose  Sand  and  Clay 
Low  200  lb  -  250  lb 

Medium  300  lb  -  350  lb 

High  400  lb  -  450  lb 

Load  Ranges  for  Compacted  Sand 
Low  300  lb  -  350  lb 

Medium  450  lb  -  500  lb 

High  700  lb  -  750  lb 

b.  Test  Procedure 

After  preparation  of  the  soil,  cone  penetrometer  measurements  were  recorded 
along  the  path  where  the  wheel  would  pass.  The  desired  load  was  then  set  up  by- 
switching  to  the  automatic  loading  mode  while  the  wheel  was  at  the  beginning  of 
the  bin.  The  load  measured  by  the  dynamometer  was  noted  on  the  recorder;  if  the 
desired  load  was  not  obtained,  air  from  a  high  pressure  supply  was  proportioned 
more  suitably  in  the  two  air  pressure  tanks  of  the  belofram  system.  The  slip 
rate  was  established  by  presetting  the  wheel  speed  at  an  estimated  condition, 
conducting  short  pilot  tests,  and  noting  the  slip  obtained.  Repeated  trials 
could  obtain  a  close  approximation  to  the  desired  slip.  The  test  was  then  ready 
to  be  conducted. 

One  person  controlled  the  wheel  and  one  person  controlled  the  recorders. 

The  person  at  the  recorders  would  first  start  the  recorders  for  a  few  seconds 
before  the  load  was  applied  for  each  channel  to  register  zero  readings.  Then 
the  chart  recording  the  sensor  outputs  was  stopped.  To  start  the  test,  the 
carriage  was  started  and  the  turning  wheel  was  lowered  into  the  fresh  soil. 
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When  the  sensors  were  about  90°  from  bottom  center,  the  chart  recording  the 
sensor  outputs  was  started;  it  was  then  turned  off  when  all  signals  of  the  100° 
sweep  were  completed.  This  was  done  to  conserve  paper  due  to  the  high  paper 
speed  required  for  clear  and  precise  signals  every  two  degrees  of  -wheel  rev¬ 
olution.  For  each  revolution,  the  paper  recording  the  sensors  was  started  and 
stopped  until  the  end  of  the  run.  After  the  run  was  complete,  both  recorders 
were  left  running  for  a  few  seconds  to  re-establish  the  zero  readings.  Cone 
penetrometer  tests  were  taken  in  the  rut  before  the  soil  was  prepared  for  the 
next  test. 

7.  Problems  Encountered 

In  concept,  this  program  appeared  to  be  straightforward  and  similar  to  that 
conducted  by  Shamay  (Reference  10)  on  the  same  equipment  in  1971.  However,  it 
was  soon  found  that  many  aspects  of  the  equipment  were  not  well  suited  to  the 
data  detail  required  by  this  program. 

The  first  problem  centered  about  the  sensors.  Initially,  the  sensors  came 
into  direct  contact  with  the  soil.  Soon  after  tests  began,  it  was  discovered 
that  sand  entered  the  spaces  between  the  sensors  and  their  surrounding  grid, 
thus  greatly  distorting  the  tangential  readings.  To  overcome  this  difficulty, 
a  rubber  membrane  was  placed  over  the  entire  circumference  of  the  wheel. 

This  membrane  solved  the  sand  intrusion  problem  but  created  others.  There 
was  a  difficulty  in  obtaining  a  good  bond  between  the  membrane  and  the  sensors 
in  order  to  get  good  transmission  of  shearing  forces.  With  the  membrane  attached 
to  the  sensors,  the  strain  in  the  rubber  became  part  of  the  sensor  system;  hence, 
each  sensor  had  to  be  calibrated  with  the  membrane  attached.  Finally,  uncer¬ 
tainty  now  arose,  when  converting  the  recorded  sensor  force  to  pressure,  as  to 
what  proportion  of  the  area  of  the  membrane  suspended  between  the  sensor  and  the 
grid  should  be  considered  as  the  bearing  surface  of  the  sensor. 

Obtaining  stable  vertical  loads  on  the  wheel  posed  a  second  problem. 
Theoretically,  the  pressure  in  the  belofram  dictated  the  load  on  the  wheel. 
Friction  in  the  system,  however,  made  the  load  vary  over  a  range  of  about 
20  pounds  during  tests. 

Proper  control  of  wheel  speed  was  also  a  problem.  At  the  extremely  slow 
speeds  required  to  separate  the  sensor  data  every  two  degrees,  the  hydraulic 
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pump/motor  system  experienced  porting  problems  that  resulted  in  a  nonuniform 
rotational  speed.  Since  the  carriage  maintained  a  constant  forward  motion,  the 
unsteady  motion  of  the  wheel  resulted  in  a  nonuniform  slip  condition.  This 
problem  would  be  avoided  by  driving  the  wheel  through  a  large  gear  reduction 
system  so  that  motor  irregularities  will  be  less  pronounced. 

During  the  testing  of  the  wheel  in  the  loam,  it  was  noted  that  an  outline 
of  the  sensors  was  visible  in  the  loam  where  the  sensors  had  been  in  contact  with 
it  (Figure  25).  Apparently  the  loam  forces  itself  and  the  membrane  into  the 
area  around  the  sensor,  thus  interfering  with  the  tangential  motion  of  the  sensor. 
The  magnitude  of  this  influence  could  not  be  determined. 

8.  Data  Processing 

After  the  tests,  a  preliminary  analysis  was  conducted  at  the  Davidson 
Laboratory  to  simplify  the  data,  convert  it  to  digital  form,  and  perform 
preparatory  validity  checks. 

All  data  except  cone  penetrometer  readings  were  initially  recorded  in 
analog  form  on  paper  strip  charts.  The  major  effort  of  data  reduction  centered 
about  the  eight  channels  of  output  from  the  four  load  sensors  mounted  on  the 
face  of  the  wheel. 

The  first  step  in  reducing  the  data  that  was  reported  from  the  sensors 
was  to  convert  the  analog  output  of  the  Sanborn  recorder  to  digital  information. 
This  conversion  was  done  manually  by  a  Gerber  scanner  at  the  Davidson  Laboratory. 
After  manually  positioning  cross  hairs  on  the  analog  curve,  the  Gerber  scanner 
would  automatically  punch  out  a  computer  card  with  the  numbers  proportional  to 
the  magnitude  recorded  at  each  of  the  eight  channels.  The  digitizing  was  con¬ 
tinued  past  bottom  center  until  all  eight  channels  had  returned  to  zero.  In 
addition  to  the  digitized  sensor  data,  each  card  contained  the  run  number,  and 
the  span  of  angular  positions  associated  with  that  card.  One  thousand  digital 
units  were  assigned  each  channel,  which  spanned  50  mm  of  recording  paper. 

Another  part  of  the  preliminary  data  reduction  scheme  was  the  computation 
of  the  applied  loads  to  the  wheel.  Using  the  output  of  the  Gerber,  and  the 
calibration  data  acquired  each  day,  a  computer  program  transferred  the  measured 
sensor  loads  at  each  interval  and  computed  the  net  horizontal  (drawbar  pull) 
and  vertical  load  and  torque.  This  computation  consisted  of  a  numerical  inte- 
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gration  of  each  2 -degree  force  measurement,  appropriately  modified  by  the  sine 
or  cosine  of  the  angle  from  the  vertical  of  the  sensors  at  the  time  of  the  mea¬ 
surement.  Consideration  was  also  given  to  the  placement  of  the  sensors  across 
the  face  of  the  wheel  and  to  the  fact  that  only  one  side  of  the  wheel  was  instrumented. 


For  further  analysis,  a  digitized  data  file  was  prepared  for  each  run, 
punched  on  paper  tape  and  transmitted  for  use  with  the  On-line  computer  system 
at  Grumman.  A  typical  printout  of  such  a  data  file  is  shown  in  Figure  26.  The 
paper  tape  was  read  at  Grumman  and  was  stored  in  the  computer  for  ready  avail¬ 
ability.  Computer  programs  were  written  at  Grumman  to  prepare  reduced  data 
files  suitable  for  graphical  display  of  the  experimental  results  on  the  visual 
display  terminal.  The  program  allowed  display  and  visual  inspection  of  the 
experimental  results  in  the  following  forms: 

•  Distribution  of  normal  stresses  across  the  wheel  at  verious  central 
angles 

•  Distribution  of  shear  stress  across  the  wheel  at  various  central  angles 

•  Longitudinal  distribution  of  normal  stresses  measured  by  individual  gages 

•  Longitudinal  distribution  of  shear  stresses  measured  by  individual  gages 

•  Longitudinal  distribution  of  interface  friction  coefficient  measured  by 

individual  sensors 

•  Longitudinal  distribution  of  average  normal  and  shear  stresses 

•  Longitudinal  distribution  of  average  interface  friction  coefficient 

In  the  preparation  of  the  reduced  data  files  for  the  display  of  average 
stresses,  the  program  allowed  the  optional  elimination  of  one  of  the  four  sensor 
readings  from  the  computation  of  averages.  This  option  made  it  possible  to  use 
data  from  some  runs  where  readings  of  a  sensor  were  found  defective. 

Copies  of  the  graphical  presentations  on  the  visual  display  terminal  were 
made  by  a  hard  copier  and  used  for  the  purposes  of  records  and  preparation  of 
illustrations. 

B.  Performance  of  Experiments  at  Grumman 
1.  Mobility  Testing  Bin 

The  Grumman  Research  Department's  mobility  testing  bin  measures  1'  x  0.85' 
x  7.5'  in  its  present  configuration.  The  bin  is  positioned  in  a  larger  container 
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and  its  width  is  adjusted  by  changing  the  size  of  the  interior  lateral  support 
members  (see  Figure  27).  This  feature  was  incorporated  into  the  bin  design  to 
facilitate  handling  of  test  bed  materials.  The  bin  can  also  be  sloped  up  to  15 
degrees,  a  feature  unique  to  the  Grumman  bin.  When  the  bin  is  sloped,  the  wheel 
assembly  can  be  pivoted  on  the  carriage  and  locked  into  place  to  allow  the  appli¬ 
cation  of  vertical  loads.  The  carriage  that  contains  the  wheel  assembly  rests  on 
teflon  sliders  and  is  driven  by  a  recirculating  ball  drive.  The  wheel  itself  is 
driven  by  a  variable  speed  motor  with  controls  that  ensure  a  constant  torque  over 
the  range  of  vertical  loads  used  in  the  testing  program.  An  aluminum  wheel,  8 
inches  in  diameter  and  2  inches  wide,  was  used  although  the  wheel  assembly  can 
accept  a  slightly  larger  diameter  and  width. 

2.  Instrumentation 

Load  and  drawbar  pull  are  obtained  directly  by  load  cells  attached  to  the 
wheel  assembly.  Normal  and  shear  stresses  are  measured  by  a  sensor  that  consists 
of  strain  gauges  on  an  axle-mounted  cantilever  beam.  The  sensor  head  (3/8-inch 
square)  is  tangent  to  the  circumference  of  the  wheel.  It  protrudes  through  an 
approximately  l/2-inch  square  opening  in  the  face  of  the  wheel.  The  opening  is 
filled  with  a  lightweight,  flexible,  felt  cloth  to  prevent  bed  material  from 
entering  the  sensor  or  lodging  in  the  sensor  housing.  As  an  added  precaution, 
a  rubber  membrane  is  stretched  over  the  wheel  so  that  the  entire  wheel  face  is 
covered.  The  center  of  the  sensor  may  be  positioned  on  the  wheel  face  either 
at  the  center  or  5/l6-inch  from  one  of  the  edges.  Torque  is  measured  directly 
by  a  Lebow  torque  sensor  on  the  shaft  of  the  wheel  drive  motor.  On  each 
revolution  of  the  wheel,  the  position  of  the  sensor  is  recorded  as  it  passes 
through  "12  o'clock."  The  travel  length  of  the  wheel  per  revolution  is 
determined  by  measuring  directly  the  distance  between  sensor  imprints  in  the 
bed  material.  A  check  is  obtained  by  comparing  the  computed  slip  to  the  slip 
corresponding  to  the  precalibrated  settings  of  the  wheel  and  carriage  motor 
controls.  The  output  signals  of  all  sensors  and  load  cells  are  transmitted 
through  a  slip  ring  on  the  wheel  axle  to  a  six-channel  Brush  strip  chart  recorder 
(Model  260). 
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V.  RESULTS  OF  EXPERIMENTS 


A.  General. 

The  experimental  data  obtained  from  each  run  (one  revolution  of  the  wheel) 
were  examined  to  determine  whether  the  measurements  were  acceptable.  The  first 
step  in  this  examination  was  a  comparison  of  the  measured  load,  drawbar  pull, 
and  torque  values  with  those  computed  by  integrating  the  measured  interface  stresses. 
In  this  comparison,  differences  can  be  expected  not  only  due  to  the  experimental 
inaccuracies,  but  also  because  of  the  approximations  in  the  integration  procedure. 

The  stresses  measured  by  the  individual  gauges  are  averaged  on  the  assumption 
that  the  stress  distribution  between  gauges  is  linear,  resulting  in  an  inherent 
inaccuracy  in  the  computation  whenever  this  is  not  the  case.  In  comparing  the 
measured  and  computed  values,  the  experimental  data  were  accepted  when  the  load 
and  torque  values  agreed  reasonably  well.  With  respect  to  the  computed  and 
measured  drawbar  pull  values,  the  criteria  for  accepting  the  test  results  were 
liberal,  mainly  because  the  measured  drawbar  pull  values  reflected  an  average 
value  over  a  full  revolution  of  the  wheel,  while  the  measured  stress  values  were 
valid  for  the  short  period  while  the  instrumented  portion  was  in  contact  with 
the  soil.  Due  to  the  problem  of  uneven  motion  of  the  wheel,  mentioned  in  Section 
IV.  A. 7,  the  difference  in  measured  and  computed  drawbar  pull  values  did  not 
necessarily  indicate  inconsistency.  In  comparison  to  theoretical  predictions, 
the  drawbar  pull  computed  from  interface  stress  measurements  was  accepted  as 
representative  of  that  particular  position  of  the  wheel. 

Another  examination  of  the  experimental  data  consisted  of  the  inspection 
of  the  transverse  distribution  of  measured  stresses.  In  some  cases,  stresses 
measured  by  an  individual  gauge  were  found  inconsistent  with  those  measured  by 
the  neighboring  gauges.  Such  defective  measurements  may  have  been  caused  by 
the  malfunction  of  electronic  circuitry,  but  may  also  have  been  caused  by  local 
soil  conditions.  Although  the  test  sand  contained  only  an  insignificant  amount 
of  sizes  greater  than  l/4-inch,  occasionally  such  a  particle  could  have  con¬ 
ceivably  been  encountered  by  a  sensor,  causing  erratic  signals.  When  such  a 
situation  was  detected,  the  defective  gauge  readings  were  eliminated  from  the 
averaging  process. 


36 


The  above  examination  of  experimental  data  resulted  in  detecting  equipment 
malfunctions  and  improper  operational  procedures,  leading  to  subsequent 
improvements  in  the  performance  of  the  experiments.  Those  experiments  where 
the  validity  of  data  was  questionable  were  eliminated  from  the  comparison  of 
experiments  to  theory  and  are  not  reported  here.  The  results  of  valid  ex¬ 
periments  are  reported  below,  grouped  according  to  the  type  of  soil  in  ■which 
they  were  performed. 


B.  Experiments  Performed  in  Loose  Sand 

The  results  of  experiments  performed  in  loose  sand  are  shown  in  Table  1, 
together  with  the  wheel  performance  data  predicted  by  the  theory.  The  measured 
distribution  of  interface  normal  and  shear  stresses  are  shown  in  Figures  27 
through  44  as  indicated  in  the  tabulation.  The  stress  distributions  for  several 
runs  that  were  performed  in  the  same  carriage  pass  are  shown  in  the  same 
illustration.  These  are: 


Fig.  30 
Fig.  31 
Fig.  36 
Fig.  37 
Fig.  38 
Fig.  44 


for  Runs  Nos.  77-78 
for  Runs  Nos.  80-81 
for  Runs  Nos.  87-88 
for  Runs  Nos.  89-91 
for  Runs  Nos.  9^-94 
for  Runs  Nos.  121-123. 


These  illustrations  show  that  the  experiments  yielded  repeatable  results 
within  the  accuracy  of  the  equipment  used  and  the  limitations  of  preparing  a 
uniform  soil  bed.  It  is  interesting  to  note  that  stress  distributions  that 
could  be  called  identical  on  visual  inspection  yield  sometimes  significant 
differences  in  drawbar  pull. 


The  measured  coefficient  of  interface  friction  was  also  evaluated  for 
most  of  the  runs.  Results,  where  available,  are  shown  in  Figs.  45  to  55.  The 
scatter  of  the  points  close  to  the  entry  and  rear  angle  is  the  result  of  both 
the  shear  and  normal  stresses  in  this  region  being  very  low.  The  inherent  in¬ 
accuracies  in  the  measurements  become  magnified  when  the  interface  friction 
coefficient  is  computed  as  the  ratio  of  shear  stress  to  normal  stress.  Further 
discussion  of  the  interface  friction  is  presented  in  Section  VI. 
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Table  1. 


RESULTS  OF.  TESTS  PERFORMED  IN  LOOSE  SAND 


Run 

Load 

Drawbar  Pull 

Torque 

Slip 

Sinkage 

Distr.  of  Avg. 

Distr.  of  Interf. 

Cone  Penetrometer 

No. 

T 

(Lbs) 

(Lbs) 

(Ft -Lb) 

a) 

(in.) 

Stresses  Fig.  No. 

Fr.  Coeff.  Fig.  No. 

Tests  Fig.  No. 

69 

M 

402 

6 

157 

33 

3.2 

45 

P 

386 

0 

107 

13 

2.5 

28 

70 

M 

446 

26 

193 

32 

3.6 

46 

P 

429 

13 

152 

20 

2.9 

71 

M 

P 

417 

386 

6 

1 

156 

108 

17 

13 

4.0 

2.5 

29 

47 

>  56 

77 

M 

519 

-1 

199 

25 

3.5 

P 

469 

4 

125 

12 

2.9 

30 

78 

M 

522 

-5 

193 

26 

3.5 

P 

474 

1 

123 

12 

2.9 

/ 

SO 

M 

4o6 

2 

146 

29 

3.6 

\ 

P 

386 

-1 

105 

12 

2.5 

31 

48 

81 

M 

379 

10 

143 

29 

3.6 

P 

385 

-1 

105 

15 

2.5 

82 

M 

P 

353 

386 

-8 

-4 

82 

90 

22 

10 

3.2 

2.3 

32 

- 

83 

M 

457 

34 

122 

14 

3.0 

33 

P 

482 

23 

199 

25 

3.5 

84 

M 

327 

-20 

50 

4 

3.5 

34 

P 

299 

-11 

45 

4 

1.6 

85 

M 

P 

319 

309 

-13 

-5 

59 

73 

4 

8 

3.0 

2.0 

35 

- 

87 

M 

P 

375 

378 

'•9 

4 

125 

122 

30 

18 

4.2 

2.6 

36 

49 

89 

M 

409 

-3 

116 

20 

4.4 

\ 

50 

P 

383 

1 

115 

15 

2.5 

f 

90 

M 

P 

44l 

455 

-2 

5 

137 

137 

21 

14 

4.4 

2.8 

\  37 

51 

91 

M 

453 

-16 

i4o 

24 

4.4 

\ 

52 

P 

465 

-3 

115 

9 

2.9 

1' 

92 

M 

413 

19 

139 

15 

3.6 

53 

P 

373 

6 

127 

18 

2.6 

93 

M 

P 

411 

380 

7 

4 

123 

120 

16 

16 

3.6  1 

2.6  | 

\  38 

54 

94 

M 

428 

13 

135 

13 

3.6 

[ 

- 

P 

444 

9 

146 

16 

2.4 

1 

102 

M 

P 

444 

440 

25 

12 

145 

153 

:  34 

1  18 

3.0 

2.9 

39 

- 

>  57 

105 

M 

532 

2 

164 

22 

3.5 

40 

P 

480 

8 

148 

14 

3.2 

117 

M 

229 

17 

68 

10 

1.4 

4i 

P 

217 

12 

70 

16 

1.8  ■ 

H9 

M 

299 

37 

113 

9 

2.0 

42 

55 

P 

316 

34 

145 

32 

2.6 

120 

M 

286 

0 

67 

5 

1.8 

\  .3 

P 

276 

6 

57 

6 

7.6 

121 

M 

250 

30 

96 

13 

1.8 

/ 

P 

248 

30 

112 

28 

2.4 

( 

122 

M 

254 

19 

83 

22 

1.8 

(  44 

P 

267 

17 

86 

16 

2.0 

123 

M 

247 

ll 

75 

21 

1.8 

I 

P 

227 

10 

68 

14 

1.7 

1 

V 

t  M  -  Measured 

P  =  Predicted  (these  values  may  differ  marginally  from  those  predicted  by  the  delivered  program 
because  of  changes  in  the  assumption  of  initial  values). 
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Predictions  shown  in  Table  1  are  based  on  cp  =  38°  friction  angle  in  the 

front  field.  In  the  rear  field,  5  percent  increase  of  the  friction  angle  was 

assumed  for  lower  loads  (up  to  300  pounds)  and  10  percent  for  high  loads.  Slip 

predictions  are  based  on  the  constants  i  =  0.07  and  K  =  0.35.  Results  of  cone 

o 

penetration  tests  are  shown  in  Figs.  56  and  57.  Parameters  in  the  plate  sinkage 
equation, 

i  n 

p  =  p  +  kz 
*o 

were  found  as  follows  (dimensions  in  inches  and  pounds): 

Po  =1.8  k  =  4.0  n  =  1.28 
C.  Experiments  Perforated  in  Dense  Sand 

The  results  of  experiments  performed  in  dense  sand  are  shown  in  Table  2, 
together  with  the  wheel  performance  data  predicted  by  the  theory.  Figures  58 
through  69  show  the  distribution  of  average  interface  stresses  obtained  in 
various  runs  for  the  loading  conditions  indicated  in  Table  2.  Figure  63  shows 
interface  stress  distributions  for  Runs  111  and  112,  which  were  obtained  in  one 
carriage  pass.  It  can  be  seen  that  the  measured  stresses  were  reasonably  well 
reproduced  in  the  two  runs,  indicating  the  validity  of  the  measurements.  The 
interface  friction  coefficients  for  these  tests  are  shown  in  Figs.  69  through 
78.  While  the  interface  friction  that  developed  along  the  soil  wheel  interface 
in  loose  sand  was  reasonably  uniform  in  most  of  the  runs,  in  dense  sand  it 
decreased  from  the  entry  and  rear  angles  toward  the  separation  angle. 

The  results  of  cone  penetrometer  tests  performed  both  before  and  after  runs 
are  shown  in  Figs.  79  through  83.  It  is  interesting  to  note  that  cone  indices 
measured  in  the  rut  were  generally  lower  than  those  measured  in  the  compacted 
soil  bed  before  the  test.  This  finding  is  consistent  with  soil  mechanics 
theory  that  associates  volume  changes  with  the  development  of  the  shear  strength 
of  granular  materials.  While  loosening  of  dense  granular  materials  on  shearing 
have  been  observed  in  many  triaxial  tests,  it  is  the  first  time  that  such  loosen¬ 
ing  has  been  observed  in  connection  with  wheel- soil  interaction.  The  wheel 
performance  predictions  given  in  Table  1  are  based  on  a  friction  angle  of 
0  =  44  degrees  in  both  the  front  and  rear  field.  Since  loosening  of  the  material 


39 


Table  2. 

.RESULTS  OF  TESTS  PERFORMED  IN  DENSE  SAND 


Run 

■ 

Load 

Drawbar  Pull 

Torque 

Slip 

Sink  age 

Dietr.  of  Avg. 

Distr.  of  Interf. 

Cone  Penetrometer 

Mo. 

fl 

(Lbs) 

(Lbs) 

( Ft- Lb) 

w 

(in.) 

Stresses  Fig.  No. 

Fr.  Coeff.  Fig.  No. 

Tests  Fig.  No. 

98 

M 

426 

138 

29 

1.5 

58 

P 

424 

158 

20 

1.5 

107 

M 

P 

486 

503 

57 

45 

158 

200 

26 

22 

1.9 

1.9 

59 

69 

108 

67 

62 

177 

240 

26 

46 

1.0 

2.2 

60 

70 

109 

56l 

530 

51 

44 

127 

205 

13 

20 

0.9 

1.8 

61 

71 

57^ 

58 

189 

no 

M 

P 

1*81 

1*67 

10 

22 

118 

116 

24 

9 

1.3 

1.2 

62 

72 

111 

M 

535 

15 

100 

30 

1.0 

P 

520 

24 

107 

6 

1.1 

73 

112 

M 

552 

4 

105 

15 

1.1 

63 

P 

567 

21 

106 

6 

1.0 

74 

\ 

113 

M 

535 

25 

124 

19 

1.1 

\ 

P 

541 

40 

188 

17 

!.7 

64 

75 

I 

n4 

M 

576 

127 

24 

1.3 

/ 

P 

542 

78 

3 

0.8 

65 

76  | 

115 

M 

513 

174 

31 

1.4 

l 

P 

527 

199 

20 

1.8 

66 

71 

/  81 

124 

M 

336 

83 

12 

1.0 

1 

P 

308 

146 

32 

1.8 

67 

78 

1 

126 

M 

550 

79 

n 

0.8 

68 

j 

P 

544 

14 

77 

3 

0.8 

i 

1 

$ 


T  M  «  Measured 

P  =  Predicted  (these  values  may  differ  marginally  from  those  predicted  by  the  delivered  program 
because  of  changes  in  the  assumption  of  initial  values). 


2+0 


is  the  end  product  of  the  shearing  process,  a  decreased  friction  angle  would  have 
manifested  itself  only  on  a  new  application  of  load,  as  in  tandem  wheel  arrange¬ 
ments  of  multiple  pass  situation.  The  slip  predictions  are  based  on  j  =  0.04 
and  K  =  0.35.  Plate  sinkage  test  parameters  for  the  compacted  sand  were  found 
as  follows: 

p  =  3.0  k  =  6.9  n  =  2 

*o 

In  some  instances,  the  predicted  torque  values  differ  appreciably  from  the 
measured  ones.  The  main  cause  of  these  discrepancies  is  that  in  the  tests  per¬ 
formed  in  dense  sand  the  distribution  of  interface  friction  was  far  from  the 
uniform  one  assumed  in  the  computations.  Further  discussion  on  the  effect  of 
non  uniform  distribution  of  interface  friction  on  wheel  performance  is  given  in 
Section  VI. 

D.  Experiments  Performed  in  Loam 

The  results  of  experiments  performed  in  loam  are  summarily  presented  in 
Table  3.  The  measured  average  noimal  and  shear  stresses  are  presented  in 
Figs.  82  through  89.  Figures  82  through  87  show  interface  stresses  measured  in 
more  than  one  run  in  the  same  pass  of  the  carriage.  The  measurements  indicate 
a  reasonable  repeatability  of  tests  run  under  the  same  conditions.  One  in¬ 
accuracy  that  occurred  in  almost  all  of  the  tests  performed  in  loam  is  a  minor 
negative  shear  stress  in  the  neighborhood  of  the  entry  angle.  This  could  have 
been  caused  either  by  some  stress  in  the  rubber  membrane  or  by  soil  intruding 
in  the  clearance  between  the  sensor  and  wheel  face.  Even  though  these  negative 
shear  stresses  were  obviously  erroneous,  their  magnitude  was  insignificant  and 
the  error  caused  by  this  inaccuracy  negligible. 

Another  interesting  feature  of  these  tests  is  the  increase  of  developed 
shear  stresses  toward  the  rear.  In  some  instances,  the  shear  stresses  in  the 
rear  were  as  high  as  the  normal  stresses  (e.g.,  Fig.  84),  indicating  a  condition 
that  is  difficult  to  explain  by  the  Mohr-Coulomb  concept  of  shear  strength. 
Examination  of  these  experimental  data  did  not  reveal  major  inconsistencies 
and,  therefore,  they  are  believed  to  be  at  least  approximately  correct.  Inter¬ 
face  stress  measurements  by  others  (Ref.  ll)  show  similar  magnitudes  of  the 
shear  stress  in  cohesive  soils,  supporting  thereby  the  validity  of  the  experi¬ 
mental  results  obtained  at  Stevens.  The  study  of  these  measurements  and  the 


Table  3 


RESULTS  OF  TESTS  PERFORMED  IN  LOAM 


Run 

Load 

Drawbar  Pull 

Torque 

Slip 

Sinkage 

Distr.  of  Avg. 

Distr.  of  Interf. 

Cone  Penetrometer 

No. 

T 

(Lbs) 

(Lbs) 

(Lb-Ft) 

(*) 

(in.) 

Stresses  Fig.  No. 

Fr.  Coeff.  Fig.  No. 

Tests  Fig.  No. 

129 

M 

336 

120 

203 

16 

1.0 

\ 

P 

336 

122 

176 

10 

1.0 

1 

) 

130 

M 

P 

322 

322 

111 

102 

185 

167 

22 

10 

1.0 

1.1 

)  82 

/  90 

131 

M 

310 

130 

199 

13 

0.9 

1 

P 

310 

109 

142 

0.6 

/ 

/ 

132 

M 

244 

24 

72 

0.9 

P 

244 

28 

28 

HE 

0.2 

91 

133 

M 

P 

244 

244 

14 

28 

59 

29 

p| 

0.9 

0.2 

83 

92 

?  97 

134 

M 

303 

100 

166 

Hill 

\ 

P 

303 

116 

92 

6 

0.4 

) 

93 

135 

M 

340 

111 

185 

7 

1.0 

V  8L 

P 

340 

103 

168 

9 

1.1 

9!* 

136 

M 

336 

114 

192 

10 

1.0 

P 

336 

112 

177 

10 

1.1 

/ 

95 

137 

M 

401 

66 

157 

0? 

1.2 

\ 

P 

401 

69 

106 

2 

0.5 

) 

- 

138 

M 

4o4 

39 

125 

0? 

wgm 

( 

P 

4o4 

39 

64 

(  85 

139 

M 

370 

34 

111 

H 

P 

370 

36 

55  • 

/ 

- 

i4o 

M 

490 

33 

l4l? 

wBM 

\ 

P 

490 

37 

70 

-1 

0.6 

) 

- 

l4l 

M 

527 

25 

149? 

-4? 

1.5 

86 

P 

508 

11 

64 

-2 

1.7 

>  98 

142 

M 

519 

0 

118? 

-8? 

1.5 

\ 

P 

479 

3 

40 

-4 

0.5 

/ 

- 

143 

M 

44l 

158 

257 

10 

1.5 

) 

P 

44l 

148 

229 

10 

1.5 

\ 

(  87 

- 

144 

M 

461 

l4o 

264 

7 

1.5 

P 

461 

128 

226 

6 

1.6 

) 

- 

145 

M 

454 

20 

121 

0? 

1.3 

P 

493 

8 

57 

-2 

0.6 

88 

- 

146 

M 

531 

107 

226 

1? 

i.3  i 

P 

531 

101 

153 

3 

0.8 

89 

96 

«• 

t  M  =  Measured 

P  «  Predicted  (these  values  may  differ  marginally  from  those  predicted  by  the  delivered  program 
because  of  changes  in  the  assumption  of  initial  values). 
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soil  properties  lead  to  the  tentative  conclusion  that  the  expansion  of  the 
partially  saturated  loam,  that  wheel  action  allows  in  this  zone,  results  in 
special  strength  properties.  Further  research  is  needed  in  this  area  to  confirm 
this  tentative  conclusion. 

This  interface  friction  coefficients  (tan  6  as  defined  in  Fig.  l)  is  shown 
for  some  runs  in  Figs.  90  through  96.  The  development  of  interface  friction 
is  different  from  that  observed  either  in  the  loose  or  dense  sand;  the  interface 
friction  appears  to  be  generally  increasing  from  the  entry  angle  towards  the  rear. 

The  predicted  values  are  based  on  a  cohesion  of  220  lbs/sq  ft  and  a 
friction  angle  of  0  =  18  degrees  as  found  by  the  triaxial  tests.  The  pre¬ 
dicted  torque  values  are  generally  lower  than  the  measured  ones  for  reasons 
explained  in  detail  in  Section  VI.  The  slip  predictions  are  based  on  the 
following  parameters:  =  0.15,  K  =  0.15.  Plate  sinkage  parameters  were 

found  to  be  p  =  4,  k  =  6.6  and  n  =  0.62. 
o 

E.  Experiments  Perforated  in  Medium  Dense  Jones  Beach  Sand 

The  results  of  experiments  perforated  at  the  Grumman  Soils  Research 
Laboratory  in  medium  dense  Jones  Beach  sand  (Long  Island,  New  York)  are  shown, 
together  with  predicted  wheel  performance  data,  in  Table  4.  The  8- inch- 
diameter,  2-inch-wide  wheel  could  accommodate  only  one  sensor  at  a  time  across 
its  width.  Therefore,  to  determine  the  stress  distribution  across  the  wheel, 
one  series  of  tests  was  performed  with  the  sensor  centerline  at  the  center  of 
the  wheel  face  and  another  series  with  its  offset  3A  inch  from  the  center. 

The  equipment  was  calibrated  prior  to  each  series  of  tests.  The  average  cone 
index  of  the  material  before  passage  of  the  wheel  was  approximately  nine  psi 
whereas  the  average  cone  index  measured  in  the  rut  after  passage  of  the  wheel  was 
approximately  15  psi.  To  account  for  this  change  of  strength,  the  predictions 
shown  in  Table  4  are  based  on  0  =  36  degrees  in  the  front  field  and  0  =  4l  degrees 
in  the  rear. 

Figures  99  and  100  show  measured  and  computed  stress  distributions  beneath 
the  wheel  for  the  level  and  sloped  bin,  respectively.  In  both  illustrations,  the 
data  points  refer  to  measured  interface  normal  and  shear  stresses,  while  the  solid 
and  dashed  curves  correspond,  respectively,  to  the  normal  and  shear  stress  distri¬ 
bution  computed  by  the  theory.  In  all  cases,  agreement  between  measured  and  com¬ 
puted  performance  parameters  is  seen  to  be  acceptable. 
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Table  4 


RESULTS  OF  TEST  PERFORMED  AT  GRUMMAN  IN  MEDIUM  DENSE  JONES  BEACH  SAND 


Run  No.* 

■ 

Load 

(Lb) 

Drawbar 

(Lbs) 

Torque 

(ft-lb) 

Slip 

(*) 

Sinkage 

(in) 

Dist.  of  . 
Fig  No 

111003 

and 

in4oi 

M 

10 

0.8 

0.9 

38 

0.53 

P 

9 

0.6 

0.5 

30 

0.37 

- 

111005 

and 

111403 

M 

10 

i 

1.1 

1.0 

17 

0.49 

_ 

P 

8 

1.6 

0.8 

16 

0.30 

“ 

110102 

M 

and 

110602 

(Level) 

7.5 

1.2 

0.7 

17 

0.32 

99 

P 

7.2 

1.8 

0.8 

21 

0.24 

99 

110104 

and 

M 

10.7 

0.8 

0.8 

17 

0.68 

100 

110604 

(Slope) 

8.9 

l.l 

1.0 

21 

0.53 

100 

T  M  =  Measured 
P  =  Predicted 


Two  run  numbers  are  given  for  each  entry,  one  corresponding 
to  the  sensor  position  at  the  middle  of  the  wheel  face  and 
the  other  for  the  sensor  at  the  edge. 
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VI.  EVALUATION  OF  TEST  RESULTS  AM)  PREDICTIONS 


A.  Validity  of  Basic  Concepts 

The  experimental  program  results  generally  confirm  the  validity  of  the 
basic  concept  of  soil-wheel  interaction  that  the  interface  stresses  are  governed 
by  the  plastic  state  of  stresses  in  the  soil.  Sample  comparisons  of  interface 
stresses  measured  and  computed  on  this  basic  concept  are  shown  in  Figs.  39>  40, 
58  and  86  for  various  conditions.  These  illustrations  show  good  agreement 
between  experiments  and  theory.  A  detailed  examination  of  all  test  results 
showed  that  there  were  three  areas  where  evaluation  of  the  experimental  results 
has  a  bearing  on  the  basic  concept.  These  are  discussed  below. 

•  The  predicted  stress  distribution  curves  from  the  front  and  rear 
field  form  a  cusp  at  the  angle  of  separation.  The  stress  distri¬ 
bution  curves  obtained  in  the  experiments  show  a  rather  smooth 
transition  at  this  point.  Whether  this  smoothness  is  due  to  some 
inertial  lag  in  the  instrumentation  or  some  adjustment  in  the  soil 
could  not  be  ascertained.  Since  the  observed  discrepancies  between 
predicted  and  measured  values  are  limited  to  a  small  area,  the 
question  is  mainly  of  academic  interest  and  does  not  affect  pre¬ 
diction  results  appreciably. 

•  In  cohesive  soils,  the  theory  predicts  an  instantaneous  rise  of  the 
interface  stresses  at  the  entry  and  rear  angles.  The  observed  stress 
rise,  though  rapid,  is  not  instantaneous.  This  discrepancy  between 
theory  and  experiments  is  the  consequence  of  using  the  Mohr-Coulomb 
yield  criterion  in  plasticity  theory  irrespective  of  the  volumetric 
strain  that  is  associated  with  the  development  of  yield  strength  of 
soil.  In  the  prediction  method,  a  stress  distribution  that  features 
instantaneous  rise  results  in  a  smaller  entry  angle  and,  consequently, 
a  lower  torque  value  than  in  reality. 

If  soil  strength  properties  are  determined  by  triaxial  tests,  the 
vertical  strain  associated  with  yield  strength  is  available.  In 
soil-wheel  interaction,  the  direction  of  major  principal  stress  in 
the  front  field  is  close  to  the  vertical,  just  as  it  is  in  the  tri¬ 
axial  test.  Thus,  an  analogous  situation  exists  from  which  approxi¬ 
mate  relationships  could  be  developed  for  the  strength  properties 
applicable  at  central  angles  close  to  the  entry  angle  where  soil  strain 


is  low.  Further  research  in  this  area  is  recommended. 

•  One  basic  assemption  of  the  theory  is  that  the  stresses  in  planes 
parallel  to  the  plane  of  motion  are  the  same,  i.e.,  the  problem 
is  two  dimensional.  An  evaluation  of  the  experimental  results  was 
made  to  see  how  well  the  experiments  conform  with  this  basic  assump¬ 
tion.  The  results  of  this  evaluation  are  summarized  below. 

In  the  tests  performed  in  loose  sand,  the  distribution  of  stresses 
across  the  wheel  was  reasonably  uniform.  A  typical  example  of  such 
distribution  is  shown  in  Figs.  101  and  102. 

In  the  tests  performed  in  dense  sand,  the  distribution  of  stresses 
across  the  wheel  was  found  to  conform  with  the  hypothesis  that  there 
is  a  limiting  transverse  distribution  governed  by  potential  lateral 
failure.  In  granular  soils,  this  limit  is  approximately  linear 
(Fig.  103)  and  varies  with  the  depth  of  the  cross  section.  Stresses 
computed  from  the  conditions  in  the  plan  of  motion  may  not  exceed 
this  limit  set  by  the  transverse  conditions.  Figures  104  and  105 
show  an  example  where  this  limiting  condition  governs  the  stresses 
across  the  full  width  of  the  wheel,  while  Figures  106  and  107  exemplify 
a  case  where  the  transverse  limiting  conditions  govern  in  that  portion 
of  the  wheel  that  is  close  to  the  side  and  longitudinal  conditions 
govern  the  stresses  in  the  center  portion. 

In  the  tests  performed  in  loam,  the  transverse  distribution  of  stresses 
was  found  to  be  influenced  by  the  condition  of  the  uniform  vertical 
displacement  imposed  on  the  soil  by  the  rigidity  of  the  wheel.  These 
conditions  require  that  the  stresses  at  the  edges  be  higher  at  the 
edges  than  in  the  center.  A  typical  transverse  distribution  of  normal 
stresses  in  loam  is  shown  in  Figs.  108  and  109.  This  type  of  distribu¬ 
tion  also  indicates  that  in  the  transverse  direction  the  soil  is  far 
from  failure  state;  although  the  hypothesis  for  a  transverse  limiting 
conditions  could  hold  true  for  cohesive  as  well  as  frictional  soils, 
in  cohesive  soils  this  condition  is  rarely  critical  as  it  will  be  seen 
from  the  discussion  below. 

The  above  typical  examples  show  various  types  of  transverse  distribution 
of  stresses  and  raise  a  question  as  to  why  the  transverse  limiting  conditions 
govern  in  one  case  and  not  in  another  one.  The  answer  lies  in  understanding 
the  nature  of  the  limiting  conditions  shown  in  Fig.  103.  The  intercept,  q, , 
is  proportional  to  the  depth  of  the  cross  section.  For  the  transverse  limit- 
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ing  conditions  to  govern,  the  sinkage  of  the  wheel  must  be  low  (depth  of 
cross  sections  small)  and  the  longitudinal  stresses  must  be  relatively  high. 

Both  of  these  conditions  are  present  in  the  wheel  tests  perfoimed  in  dense 
sand  and  are  absent  in  the  test  perfoimed  in  loose  sand. 

In  cohesive  soils,  the  magnitude  of  q^.  intercept  depends  not  only  on 
the  depth  of  the  cross  section,  but  also  on  the  value  of  cohesion.  Even  a 
small  amount  of  cohesion  is  sufficient  to  result  in  such  q_^_  values  that  the 
transverse  limiting  condition  ceases  to  be  critical. 

It  is  possible  to  foimulate  the  above  qualitative  statements  mathemati¬ 
cally  and  to  improve  thereby  the  accuracy  of  the  prediction  method.  A  further 
advantage  of  the  mathematical  formulation  would  be  in  the  analysis  of  multiple 
wheel  perfoimance.  Adjacent  wheels  influence  the  transverse  limiting  conditions 
by  hindering  lateral  failure;  their  effect  would,  however,  be  negligible  when 
transverse  limiting  conditions  are  not  critical. 

B.  Validity  of  Tentative  Assumptions 

In  the  application  of  the  basic  concept  of  soil-wheel  interaction  to 
the  problem  of  wheel  perfoimance  calculations,  it  was  necessary  to  make  certain 
assumptions,  discussed  in  Section  II,  to  define  the  problem  completely.  These 
assumptions  were  made  on  the  basis  of  experimental  information  available  at  the 
time  of  the  development  of  the  program.  With  more  information  available  from 
the  validation  test  series  performed  at  Stevens  Institute,  it  is  proper 
to  reexamine  these  assumptions.  The  results  of  this  reexamination  are  summar¬ 
ized  below. 

For  the  definition  of  the  boundary  condition  at  the  interface,  the  inter¬ 
face  friction  angle  6  was  introduced.  In  the  computations,  this  angle  was 
assumed  to  be  constant.  This  assumption,  however,  is  not  essential  to  the 
basic  concept  that  requires  only  that  the  angle  6  be  defined  but  not  necessar¬ 
ily  constant  along  the  interface.  For  the  evaluation  of  the  interface  friction 
that  developed  along  the  interface  in  the  experiments,  a  computer  program  was 
written  that  calculates  the  coefficient  of  interface  friction  (tan  6)  from 


the  raw  data  file  and  prepares  a  data  file  suitable  for  viewing  on  the  visual 
display  terminal.  The  computed  values  are  averages  over  the  width  of  the 
wheel  and  are  shown,  where  available,  in  the  illustrations  referenced  in 
Tables  1  through  3.  The  following  general  conclusions  may  be  drawn  from 
these  illustrations: 

In  loose  sand,  with  some  exceptions,  the  assumption  of  a  constant  inter¬ 
face  friction  coefficient  appears  to  be  a  reasonable  approximation  (see  Figs. 
45  through  55).  The  exceptions  are  where  the  interface  friction  coefficient 
decreases  from  the  edges  toward  the  center,  a  distribution  that  was  found 
typical  of  the  tests  performed  in  dense  sand  (see  Figs.  69  through  78). 

In  the  loam,  the  distribution  of  the  interface  friction  coefficient  along 
the  interafce  was  found  to  be  different  from  either  that  typical  of  loose 
sand  or  that  typical  of  dense  sand.  The  typical  feature  of  the  distributions 
of  the  interface  friction  coefficient  in  loam  shown  in  Figs.  90  through  98  is 
that  the  highest  value  occurs  at  =  0  degree,  the  bottommost  part  of  the 
wheel.  An  explanation  of  this  feature  could  be  that  the  developnent  of  inter¬ 
face  friction  in  the  highly  compressible  loam  is  associated  with  volumetric 
straining  of  the  soil,  which  is  obviously  the  greatest  at  the  bottom  of  the 
wheel . 

In  anticipation  of  a  nonuniform  distribution  of  the  interface  friction 
coefficient,  provisions  were  made  in  the  computer  program  that  allow  a  linear 
variation  of  the  interface  friction  coefficient  along  the  interface.  However, 
the  new  experimental  evidence  obtained  in  the  validation  test  series  is  not 
sufficient  to  formulate  a  relationship  among  the  various  types  of  the  distri¬ 
butions  of  interface  friction  coefficients  and  soil  properties  and  wheel 
loadings.  For  this  reason,  it  was  not  possible  to  use  that  capacity  of  the 
program  that  allows  a  linear  variation  of  6  along  the  interface.  Further 
theoretical  and  experimental  research  is  needed  to  clarify  the  relationships 
that  govern  the  development  of  interface  friction.  This  is  all  the  more  im¬ 
portant  since  the  traction  developed  by  wheels  is  directly  related  to  the 
developnent  of  interface  friction. 

Another  reason  to  do  further  research  in  this  area  is  the  connection 
between  the  angle  6  and  the  angle  of  separation,  ot .  In  Section  II,  a 
relationship  between  6  and  was  established  that  was  incorporated  in  the 
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computer  program.  In  the  development  of  this  relationship,  6  was  assumed  to 
be  constant.  If  6  is  variable  along  the  interface,  may  become  indefinite. 
Another  possibility  that  would  have  to  be  investigated  is  that  the  variation 
of  the  interface  friction  reflects  the  mutual  adjustment  of  the  angle  of 
separation  and  interface  friction  so  that  the  requirements  set  forth  in 
Eq.  (l6)  may  be  met.  A  review  of  predictions  indicates  that  in  most  cases 
where  prediction  accuracy  was  not  good,  the  maximum  normal  stress  occurred 
at  an  angle  that  deviates  from  the  hypothesized  separation  angle.  Thus,  a 
significant  improvement  in  prediction  accuracy  could  be  obtained  if  the 
variation  of  the  friction  along  the  interface  and  the  angle  of  separation 
associated  with  this  variation  could  be  more  accurately  introduced  in  the 
program. 

The  development  of  interface  friction  is  also  associated  with  slip  and 
Eq.  (13)  was  used  in  the  program  to  compute  slip  from  the  value.  Obviously, 
Eq.  (13)  is  not  defined  for  a  variable  6  and  this  may  be  one  reason  for  the 
poor  slip  predictions  shown  in  the  tabulations.  Unfortunately,  the  measured 
slip  values  are  also  somewhat  uncertain  because  of  the  uneven  rotation  of  the 
wheel  discussed  in  Section  IV.  A.  1.  For  this  reason,  it  is  difficult  to 
evaluate  the  validity  of  Eq.  (13)  or  to  draw  conclusions  about  the  influence 
of  a  variable  6  on  slip.  An  interesting  concept  that  would  evolve  from  the 
study  of  the  variation  of  6  and  slip  is  shown  schematically  in  Fig.  110.  It 
was  found  that,  at  lease  in  the  front  field,  lines  drawn  perpendicular  to  the 
direction  of  the  major  prinicpal  stress  intersected  the  vertical  of  the  wheel 
axle  within  very  narrow  limits  presumably  centering  around  the  instantaneous 
center  of  reaction.  Since  the  direction  of  principal  stresses  coincides  with 
that  of  principal  strains  in  isotropic  soils,  it  is  reasonable  to  assume  that 
the  displacement  velocity  vector  would  be  directed  the  same  way  as  the  princi¬ 
pal  strain  vector  in  compressible  media.  Could  this  be  proven  and  formulated 
mathematically,  a  very  important  breakthrough  in  the  somewhat  nebulous  re¬ 
lationship  between  slip  and  mobilization  of  interface  friction  could  be  accom¬ 
plished. 

To  study  possibilities  of  improving  slip  predictions,  an  analysis  of 
measured  and  predicted  slip  values  was  made.  Results  of  this  analysis  are 
shown  in  Figs.  Ill  and  112,  which  show  measured  and  predicted  slip  values 
plotted  against  the  pull  coefficient  for  the  tests  performed  in  loose  sand. 


For  a  given  soil,  the  points  should  be  within  a  narrow  band  representing  a 
unique  relationship  for  the  experimental  scatter.  It  can  be  seen  from  Fig. 

Ill  that  the  experimental  results  do  not  collapse  in  a  narrow  band,  indica¬ 
ting  that  the  uneven  motion  of  wheel  resulted  in  inaccurate  slip  measurements. 
Because  of  the  uncertainities  in  the  measured  slip  values,  it  was  pointless 
to  try  improvements  in  the  theoretical  predictions  that  yield,  at  least  quali¬ 
tatively,  a  satisfactory  pull  coefficient  slip  relationship  (Fig.  112). 
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VII  CONCLUSIONS  AND  RECOMMENDATIONS 


The  theoretical  and  experimental  investigations  performed  under  the  contract 
conclusively  show  that  the  proposed  concept  of  soil -wheel  interaction  is  valid 
and  that  the  application  of  plasticity  theory  to  wheel-tire  interaction  problems 
is  a  valuable  tool  in  the  mathematical  formulation  of  the  problem.  The  analysis 
of  experimental  results  in  the  framework  of  the  basic  concept  provided  new  in¬ 
sight  into  the  interaction  problem  and  essential  new  information  was  gained  for 
the  theoretical  formulation  of  more  complex  interaction  problems  such  as  soil 
interaction  with  pneumatic  tires,  tandem  and  multiple  wheels,  and  multipass 
interaction  analysis. 

The  computer  program  developed  for  the  numerical  solution  of  the  interaction 
problem  yields  the  answers  within  an  acceptable  computer  time .  Predictions  by  the 
computer  program  were  generally  good.  The  accuracy  of  predictions  depended  on 
how  well  certain  assumptions,  made  in  the  development  of  the  program  on  the  basis 
of  experimental  information,  approximated  actual  conditions. 

The  areas  where  further  research  would  result  in  improved  prediction  accuracy 
or  in  expanding  the  applicability  of  the  concept  and  program  to  cases  not  covered 
in  the  present  study  are  listed  below. 

•  Theoretical  and  experimental  research  in  the  area  of  the  development 
of  the  interface  friction  and  its  relationship  with  slip.  Research 

in  this  area  would  lead  to  a  significant  improvement  in  the  prediction 
method.  Also,  theoretical  formulation  of  the  development  of  interface 
friction  for  towed  and  braked  wheels  could  be  included  in  this  research 

i 

and  used  in  the  application  of  the  present  computer  program  for  the 
prediction  of  towed  and  braked  wheel  performance.  This  would  then  be 
used  to  predict  the  performance  of  2x4,  4x6  and  other  multi-axle  vehicles. 

•  Theoretical  and  experimental  research  to  develop  a  theory  for  changes 
in  strength  properties  of  soils  due  to  the  compacting  effect  of 
wheels.  Mobility  predictions  are  based  on  virgin  ground  soil  properties, 
yet  traction  develops  primarily  in  the  rear  portion  of  the  wheel  where 
the  soil  is  already  compacted.  A  study  to  formulate  the  effect  of 
compaction  would  be  essential  not  only  for  application  in  the  present 
theory  but  also  for  use  with  any  other  predictive  method.  In  the  frame¬ 
work  of  such  a  study,  an  analysis  of  the  volumetric  strain  necessary  to 
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develop  the  Mohr- Coulomb  yield  strength  could  be  made  and  improvements 
in  the  accuracy  of  wheel  performance  predictions  could  be  achieved. 
Research  in  the  strength  properties  would  also  be  useful  for  the 
formulation  of  the  analysis  of  tandem  wheel  performance  and  for  the 
development  of  multipass  criteria. 

•  Theoretical  and  experimental  research  to  formulate  criteria  for 
lateral  failure.  Such  research  would  not  only  improve  the  prediction 
method,  but  would  also  lay  the  foundations  for  the  analysis  of  multiple 
wheel  performance. 

•  Validation  test  program  for  slopes.  All  validation  tests,  except  for 
one  at  Grumman,  were  performed  on  level  soil  beds.  It  would  be  desir¬ 
able  to  perform  a  test  series  on  slopes  with  interface  stress  measure¬ 
ments  since  no  such  information  is  available  at  present. 

•  An  extension  of  the  theory  of  rigid  wheel-soil  interaction  to  pneumatic 
tire-soil  interaction  so  that  the  effect  of  tire  deflection  on  soil 
response  could  be  taken  into  account  and  a  computer  program  for  the 
prediction  of  tire  performance  in  soft  soil  formulated. 
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Fig.  4  Effect  of  the  Interface  Friction  Angle  on  the  Geometry 
of  Rear  Slip  Line  Field  and  Associated  Normal  Stresses 
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Fig.  5  Slip  Line  Fields  and  Associated  Normal  Stresses  for 
Various  Entry  Angles 
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Fig.  7  Computation  of  Slip  Line  Field  by  Finite  Differences 


INPUT:WHEEL  PARAMETERS  R,  B,  a  e,  a  m 
SOIL  PROPERTIES  C,  $  ,  Y 
TERRAIN  PARAMETER  * 

INTERACTION  PARAMETER  8 


Fig.  8  Flow  Diagram  for  the  Computation  of  a  Single 
Slip  Line  Field  (Subroutine  1,  SLIP) 
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IN  PUT:  WHEEL  PARAMETERS  R,  B 
SOIL  PROPERTIES  C,*,  Y 
TERRAIN  PARAMETER  « 
INTERACTION  PARAMETER  8 


Fig.  9  Flow  Diagram  for  the  Computation  of  a  Matching 
Set  of  Slip  Line  Fields  and  Wheel  Performance 
Parameters 
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Fig  10  Flow  Diagram  for  the  Iterative  Procedure  for  Finding  a  and  6 
Values  for  Given  Load,  Drawbar  Pull,  and  Soil  Properties 


Fig.  11  Flow  Diagram  for  the  Computation  of  Wheel  Performance 
in  Cohesive  Soils 
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Fig.  12  Estimation  of  Strength  Properties  of  Soil  in  the 
Rear  Field 


Fig.  13  L  Shaped  Sensors  Mounted  on  the 
Circumference  of  the  Wheel 
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TEST  NO  INITIAL  DENSITY 

39  1.590IG/CU  CM)  98.4ILBS/CU  FT) 
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Fig.  17  Results  of  Triaxial  Tests  on  Loose  Sand 
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Fig.  18  Results  of  Triaxial  Tests  on  Dense  Sand 


Fig.  19  Grain  Size  Distribution  of  Test  Sand 
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Fig.  20  Principal  Stress  Difference  Versus  Strain  --  Results 
of  Triaxial  Tests  in  Loam 
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Results  of  Cone  Penetrometer  Tests  in  Sand  at 
Various  Densities 
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Fig.  26  Typical  Digitized  Raw  Data  File 
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p  *  normal  stress  Run  #69 

»  +  shear  stress  Run  #69 

o  normal  stress  Run  #70 
8  0  x  shear  stress  Run  #70 


Central  Angle  a 


Fig.  28  Run  #69  - 


70  Average  Interface  Normal  and  Shear  Stresses 
Measured  in  Loose  Sand 


Run  #69 


Load  =  402  Lbs,  Drawbar  Pull  =  6  lbs, 
Torque  =  157  ft  lbs 


Run  #70 


Load  =  446  lbs,  Drawbar  Pull  =  26  lbs, 
Torque  =  193  ft  lbs 
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Fig.  29  Run  #71  Average  Interface  Normal  and  Shear  Stresses 

Measured  in  Loose  Sand 

Load  =  417  lbs.  Drawbar  Pull  =  6  lbs,  Torque  =  156  ft  lbs 
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Interface  Stresses  (lbs/sq  ft) 
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*  normal  stress  Run  #77 
+  shear  stress  Run  #77 

o  normal  stress  Run  #78 
o  x  shear  stress  Run  #78 


'10°  Central  Angle  a  40° 


Fig.  30  Run  #77-78  Average  Interface  Normal  and  Shear  Stresses 

Measured  in  Loose  Sand 

Run  #77:  Load  =  519  lbs,  Drawbar  Pull  =  1  lb, 
Torque  =  199  ft  lbs 

Run  #78:  Load  =  522  lbs,  Drawbar  Pull  =  -5  lbs, 
Torque  =  193  ft  lbs 
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Interface  Stresses  (lbs/sq  ft) 
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+  shear  stress  Run  #80 
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,  o  normal  stress  Run  #81 

,  x  shear  stress  Run  #81 


- 1  C 0  Central  Angle  a  40° 


Fig.  31  Run  #80-81  Average  Interface  Normal  and  Shear  Stresses 

Measured  in  Loose  Sand 

Run  #80:  Load  =  406  lbs,  Drawbar  Pull  =  2  lbs, 
Torque  =  146  ft  lbs 

Run  #81:  Load  =  379  lbs.  Drawbar  Pull  =  10  lbs, 
Torque  =  143  ft  lbs 
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Fig.  33  Run  #83  Average  Interface  Normal  and  Shear  Stresses 

Measured  in  Loose  Sand 

Load  =  457  lbs,  Drawbar  Pull  =  34  lbs,  Torque  =  122  ft  lbs 
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-10°  Central  Angle  a  40° 


Fig.  34  Run  #84  Average  Interface  Normal  and  Shear  Stresses 

Measured  in  Loose  Sand 

Load  =  327  lbs,  Drawbar  Pull  =  -20  lbs,  Torque  =  50  ft  lbs 


Interface  Stresses  (lbs/sq  ft) 


GRUTIWN  RESEARCH  TIME-SHARED  GRAPHICS  TERM  INAL 


Fig.  35  Run  #85  Average  Interface  Normal  and  Shear  Stresses 

Measured  in  Loose  Sand 

Load  =  319  lbs,  Drawbar  Pull  =  -13  lbs,  Torque  =  59  ft  lbs 


Interfaces  Stresses  (lbs/sq  ft) 


I  I  * 

X  X  x 


~  normal  stress  Run  #8 
+  shear  stress  Run  #8 

x  normal  stress  Run  #8 
x  shear  stress  Run  #8 


1 1 1 +  *  t 

X  ++++XxXx  X 1 

+  „  V  X  X  ■ 

t  tyr 
X  t  x 


-10°  Central  Angle  a 

36  Run  #87-88  Average  Interface  Norm 

Measured  in  Loos  e  Sand 
Run  #87:  Load  =  375  1 
Torque  =  125 
Run  #88:  Load  =  354  1 
Torque  =  100 
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Fig.  37  Run  #89-91  Average  Interface  Normal  and  Shear  Stresses 

Measured  in  Loose  Sand 

Run  #89:  Load  =  409  lbs,  Drawbar  Pull  =  -3  lbs, 
Torque  =  116  ft  lbs 

Run  #90:  Load  =  441  lbs,  Drawbar  Pull  =  -2  lbs, 
Torque  =  137  ft  lbs 

Run  #91:  Load  =  453  lbs,  Drawbar  Pull  =  -16  lbs, 
Torque  =  140  ft  lbs 
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Fig.  41  Run  #117  Average  Interface  Normal  and  Shear  Stresses 

Measured  in  Loose  Sand 

Load  =  229  lbs,  Drawbar  Pull  =  17  lbs, 

Torque  =  68  ft  lbs 
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Interface  Stresses  (lbs/sq  ft) 


Fig.  42  Run  #119  Average  Interface  Normal  and  Shear  Stresses 

Measured  in  Loose  Sand 

Load  *  299  lbs,  Drawbar  Pull  =  37  lbs, 

Torque  =  113  ft  lbs 
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Interface  Stresses  (lbs/sq  ft) 
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Fig.  43  Run  #120  Average  Interface  Normal  and  Shear  Stresses 

Measured  in  Loose  Sand 

Load  *  286  lbs,  Drawbar  Pull  «  0  lbs. 

Torque  =  67  ft  lbs 
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Fig.  44  Run#121-123  Average  Interface  Normal  and  Shear  Stresses 

Measured  in  Loose  Sand 

Run  #121:  Load  «  250  lbs,  Drawbar  Pull  =  30  lbs, 
Torque  =  96  ft  lbs 

Run  #122:  Load  ■=  254  lbs,  Drawbar  Pull  *  19  lbs, 
Torque  *  83  ft  lbs 

Run  #123:  Load  -  247  lbs,  Drawbar  Pull  -  11  lbs, 
Torque  »  75  ft  lbs 


-10° 


Central  Angle  a 


40 


Fig.  45  Run  #69  Variation  of  Interface  Friction  Coefficient 

(tan  6)  Along  the  Interface  in  Loose  Sand 
Load  =  402  lbs,  Drawbar  Pull  =  6  lbs. 

Torque  =  157  ft  lbs 


Interface  Friction  Coefficient  (tan  6) 
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-10°  Central  Angle  a  40° 


Fig.  46  Run  #70  Variation  of  Interface  Friction  Coefficient 

(tan  6)  Along  the  Interface  in  Loose  Sand 
Load  =  446  lbs,  Drawbar  Pull  =  26  lbs, 
Torque  =  193  ft  lbs 


Interface  Friction  Coefficient  (tan  5) 
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Fig.  47  Run  #71  Variation  of  Interface  Friction  Coefficient 

(tan  5)  Along  the  Interface  in  Loose  Sand 
Load  =  417  lbs,  Drawbar  Pull  =  6  lbs, 

Torque  =  156  ft  lbs 


98 


Interface  Friction  Coefficient  (tan  6) 
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Fig.  48  Run  #80-81  Variation  of  Interface  Friction  Coefficient 

(tan  6)  Along  the  Interface  in  Loose  Sand 
Run  #80:  Load  =  406  lbs,  Drawbar  Pull  =  2  lbs, 
Torque  =  146  ft  lbs 

Run  #81:  Load  =  379  lbs,  Drawbar  Pull  =  10  lbs, 
Torque  =  143  ft  lbs 
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Interface  Friction  Coefficient  (tan  5) 
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Fig.  49  Run  #87  Variation  of  Interface  Friction  Coefficient 

(tan  6)  Along  the  Interface  in  Loose  Sand 
Load  =375  lbs,  Drawbar  Pull  =  9  lbs, 

Torque  =  125  ft  lbs 
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Interface  Friction  Coefficient  (tan  6) 


Interface  Friction  Coefficient  (tan  5) 
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-10°  Central  Angle  a  40° 


Fig.  51  Run  #90  Variation  of  Interface  Friction  Coefficient 

(tan  6)  Along  the  Interface  in  Loose  Sand 
Load  =  441  lbs,  Drawbar  Pull  =  -2  lbs, 
Torque  =  137  ft  lbs 
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Interface  Friction  Coefficient  (tan  5) 
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Fig.  52  Run  #91  Variation  of  Interface  Friction  Coefficient 

(tan  5)  Along  the  Interface  in  Loose  Sand 
Load  =  453  lbs,  Drawbar  Pull  =  -16  lbs, 
Torque  =  140  ft  lbs 
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Interface  Friction  Coefficient  (tan  6) 
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Fig.  53  Run  #92  Variation  of  Interface  Friction  Coefficient 

(tan  6)  Along  the  Interface  in  Loose  Sand 
Load  =  413  lbs,  Drawbar  Pull  =  19  lbs, 
Torque  =  139  ft  lbs 
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Interface  Friction  Coefficient  (tan  5) 
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Fig.  54  Run  #93  Variation  of  Interface  Friction  Coefficient 

(tan  6)  Along  the  Interface  in  Loose  Sand 
Load  =  411  lbs,  Drawbar  Pull  =  7  lbs 
Torque  =  123  ft  lbs 
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Interface  Friction  Coefficient  (tan  6) 


GRUMMAN  RESEARCH  TIME-SHARED  GRAPHICS  TERMINAL 


Fig.  55  Run  #119  Variation  of  Interface  Friction  Coefficient 

(tan  6)  Along  the  Interface  in  Loose  Sand 
Load  =  229  lbs,  Drawbar  Pull  =  37  lbs, 
Torque  =  68  ft  lbs 
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Fig.  56  Results  of  Cone  Penetrometer  Tests  in  Loose  Sand 
for  Runs  #69-78 


Fig.  57  Results  of  Cone  Penetrometer  Tests  in  Loose  Sand 
for  Runs  #80-123 
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Fig.  59  Run  #107  Average  Interface  Normal  and  Shear  Stresses 

Measured  in  Dense  Sand 

Load  =  486  lbs,  Drawbar  Pull  =  51  lbs, 

Torque  =  158  ft  lbs 
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Fig.  60  Run  #108  Average  Interface  Normal  and  Shear  Stresses 

Measured  in  Dense  Sand 

Load  =  487  lbs,  Drawbar  Pull  =  67  lbs, 

Torque  =  177  ft  lbs 
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Fig.  61  Run  #109  Average  Interface  Normal  and  Shear  Stresses 

Measured  in  Dense  Sand 

Load  =  561  lbs,  Drawbar  Pull  =  51  lbs, 

Torque  =  127  ft  lbs 
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Fig.  62  Run  #110  Average  Interface  Normal  and  Shear  Stresses 

Measured  in  Dense  Sand 

Load  =  481  lbs,  Drawbar  Pull  =  10  lbs, 

Torque  =  118  ft  lbs 
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Fig.  63  Run  #111-112  Average  Interface  Normal  and  Shear  Stresses 

Measured  in  Dense  Sand 

Run  #111:  Load  =  535  lbs,  Drawbar  Pull  =  15  lbs, 
Torque  =  100  ft  lbs 

Run  #112:  Load  =  552  lbs,  Drawbar  Pull  =  4  lbs, 
Torque  =  105  ft  lbs 
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Fig.  64  Run  #113  Average  Interface  Normal  and  Shear  Stresses 

Measured  in  Dense  Sand 
Load  =  535  lbs,  Drawbar  Pull  =  25  lbs, 
Torque  =  124  ft  lbs 
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Fig.  65  Run  #114  Average  Interface  Normal  and  Shear  Stres 

Measured  in  Dense  Sand 

Load  =  576  lbs.  Drawbar  Pull  =  2  lbs, 

Torque  =  127  ft  lbs 
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Fig.  66  Run  #115  Average  Interface  Normal  and  Shear  Stresse 

Measured  in  Dense  Sand 

Load  =  513  lbs,  Drawbar  Pull  =  41  lbs, 

Torque  =  174  ft  lbs 
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Fig,  67  Run  #124  Average  Interface  Normal  and  Shear  Stresses 

Measured  in  Dense  Sand 

load  *  336  lbs,  Drawbar  Pull  »  42  lbs, 

Torcjue  «=  83  ft  lbs 
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Fig.  68  Run  #126  Average  Interface  Normal  and  Shear  Stresses 

Measured  in  Dense  Sand 

Load  =  550  lbs,  Drawbar  Pull  =  1  lb, 

Torque  =  79  ft  lbs 
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Interface  Friction  Coefficient  (tan  6) 


Fig.  69  Run  #107  Variation  of  Interface  Friction  Coefficient 

(tan  6)  Along  the  Interface  in  Dense  Sand 
Load  *=  480  lbs,  Drawbar  Pull  e=  51  lbs 
Torque  =  158  ft  lbs  ? 
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Interface  Friction  Coefficient  (tan  6) 
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Fig.  72  Run  #110  Variation  of  Interface  Friction  Coefficient 

(tan  6)  Along  the  Interface  in  Dense  Sand 
Load  =481  lbs,  Drawbar  Pull  =  10  lbs, 
Torque  =  118  ft  lbs 
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Interface  Friction  Coefficient  (tan  5) 
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Fig.  73  Run  #111  Variation  of  Interface  Friction  Coefficient 

(tan  6)  Along  the  Interface  in  Dense  Sand 
Load  -  535  lbs,  Drawbar  Pull  =  15  lbs, 
Torque  «=  100  ft  lbs 
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Fig.  74  Run  #112  Variation  of  Interface  Friction  Coefficient 

(tan  6)  Along  the  Interface  in  Dense  Sand 
Load  =  576  lbs,  Drawbar  Pull  =  2  lbs, 

Torque  =  105  ft  lbs 
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Interface  Friction  Coefficient  (tan  5) 


Interface  Friction  Coefficient  (tan  6) 
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Fig.  76  Run  #114  Variation  of  Interface  Friction  Coefficient 

(tan  6)  Along  the  Interface  in  Dense  Sand 
Load  =  576  lbs,  Drawbar  Pull  =  2  lbs, 

Torque  =  127  ft  lbs 
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Fig.  77  Run  #115  Variation  of  Interface  Friction  Coefficient 

(tan  6)  Along  the  Interface  in  Dense  Sand 
Load  =  513  lbs,  Drawbar  Pull  m  41  lbs. 
Torque  *»  174  ft  lbs 


128 


Interface  Friction  Coefficients  (tan  6) 


3SUMW  RESEARCH  T  IMt-SHfiRED  GRAPHICS  iOTIINAL 


Fig.  78  Run  #126  Variation  of  Interface  Friction  Coefficient 

(tan  5)  Along  the  Interface  in  Dense  Sand 
Load  =  550  lbs,  Drawbar  Pull  =  1  lb, 

Torque  =  79  ft  lbs 
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Fig.  79  Results  of  Cone  Penetrometer  Tests  in  Dense  Sand 
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Fig.  82  Run  #129-131  Average  Interface  Normal  and  Shear  Stresses 

Measured  in  Loam 

Run  #129:  Load  =  336  lbs,  Drawbar  Pull  =  126  lbs, 
Torque  =  203  ft  lbs 

Run  #130:  Load  =  322  lbs,  Drawbar  Pull  =  111  lbs, 
Torque  =  185  ft  lbs 

Run  #131:  Load  =  310  lbs,  Drawbar  Pull  =  130  lbs, 
Torque  =  199  ft  lbs 
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o  normal  stress,  Run  #132 
x  shear  stress,  Run  #132 
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Fig.  83  Run  #132-133  Average  Interface  Normal  and  Shear  Stresses 

Measured  in  Loam 

Run  #132:  Load  »  244  lbs,  Drawbar  Pull  =  24  lbs 
Torque  «=  72  ft  lbs  * 

Run  #133:  Load  «=  244  lbs,  Drawbar  Pull  «  14  lbs 
Torque  -  59  ft  lbs 
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Fig.  84  Run  #134-136  Average  Interface  Normal  and  Shear  Stresses 

Measured  in  Loam 

Run  #134:  Load  =  303  lbs,  Drawbar  Pull  =  100  lbs, 
Torque  =  166  ft  lbs 

Run  #135:  Load  =  340  lbs,  Drawbar  Pull  =  111  lbs, 
Torque  =  185  ft  lbs 

Run  #136:  Load  =  336  lbs,  Drawbar  Pull  =  114  lbs, 
Torque  «  192  ft  lbs 
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Fig.  85  Run  #137-139  Average  Interface  Normal  and  Shear  Stresses 
v  Measured  in  Loam 

Run  #137:  Load  *  401  lbs.  Drawbar  Pull  «  66  lbs, 
Torque  =  157  ft  lbs 

Run  #138:  Load  -  404  lbs,  Drawbar  Pull  m  39  lbs, 
Torque  =  125  ft  lbs 

Run  #139:  Load  *=  370  lbs,  Drawbar  Pull  *  34  lbs, 
Torque  -  111  ft  lbs 
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Fig.  86  Run  #140-142  Measured  and  Predicted  A\r 

and  Shear  Stresses'  in  Loam^86  nterface  Normal 
Run  #140:  Load  =  490  lbs,  Drawbar  PuU  =  33 

p,  7/1/1  °rqUe  =  141  ft  lbs  '  : 

n  •  Load  -  527  lbs,  Drawbar  Pull  =  25  lbs 

Run  #142-  l0olT-  nbS 

Load  -519  lbs,  Drawbar  Pull  =  0  lb 

Torque  =  118  ft  lbs  ’ 
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Fig,  87  Run  #143-144  Average  Interface  Normal  and  Shear  Stresses 

Measured  in  Loam 

Run  #143:  Load  =441  lbs,  Drawbar  Pull  -  158  lbs, 
Torque  =  251  ft  lbs 

Run  #144:  Load  =  461  lbs,  Drawbar  Pull  =  140  lbs, 
Torque  =  226  ft  lbs 
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Fig.  88  Run  #145  Average  Interface  Normal  and  Shear  Stresses 

Measured  in  Loam 

Load  =  454  lbs,  Drawbar  Pull  =  20  lbs. 
Torque  =  121  ft  lbs 
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Fig.  90  Run  #129-131  Variation  of  Interface  Friction  Coefficient 

(tan  5)  Along  the  Interface  in  Loam 
Run  #129:  Load  =  336  lbs,  Drawbar  Pull  =  126  lbs, 
Torque  =  203  ft  lbs 

Run  #130:  Load  -  322  lbs,  Drawbar  Pull  =  111  lbs, 
Torque  =  185  ft  lbs 

Run  #131:  Load  =  310  lbs,  Drawbar  Pull  =  130  lbs, 
Torque  =  199  ft  lbs 
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Fig.  92  Run  #133  Variation  of  Interface  Friction  Coefficient 

(tan  6)  Along  the  Interface  in  Loam 
Load  =  244  lbs,  Drawbar  Pull  =  14  lbs, 
Torque  =  59  ft  lbs 
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Fig.  93  Run  #134  Variation  of  Interface  Friction  Coefficient 

(tan  6)  Along  the  Interface  in  Loam 
Load  303  lbs,  Drawbar  Pull  =  100  lbs, 
Torque  166  ft  lbs 


ibk 


Interface  Friction  Coefficient  (tan  6) 


mm  RESEARCH  Ti I1E -SHARED  GRAPHICS  TERM  NHL 


Fig.  94  Run  #135  Variation  of  Interface  Friction  Coefficient 

(tan  8)  Along  the  Interface  in  Loam 
Load  =  340  lbs,  Drawbar  Pull  =  111  lbs, 
Torque  =  185  ft  lbs 


145 


Interface  Friction  Coefficient  (tan  6) 


mm  RESEARCH  TIME ‘SHARED  GRAPHICS  TERMINAL 


-Iff  Central  Angle  a  X  40° 


Fig.  95  Run  #136  Variation  of  Interface  Friction  Coefficient 

(tan  6)  Along  the  Interface  in  Loam 
Load  ■  336  lbs,  Drawbar  Pull  -  114  lbs, 
Torque  -  192  ft  lbs 
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Fig.  96  Run  #146  Variation  of  Interface  Friction  Coefficient 

(tan  6)  Along  the  Interface  in  Loam 
Load  =  531  lbs,  Drawbar  Pull  =  107  lbs, 
Torque  =  226  ft  lbs 
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Fig.  98  Results  of  Cone  Penetrometer  Tests  in  Loam  for 
Runs  #137-146 
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Fig.  99  Run  #110102  Measured  and  Predicted  Average  Interface 

Normal  and  Shear  Stresses.  Wheel 
Diameter  ■  8  in.,  Width  =  2  in.  -  Jones 
Beach  Sand  -  Level  Surface 
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Fig.  100  Run  #110602  Measured  and  Predicted  Average  Interface 

Normal  and  Shear  Stresses.  Wheel 
Diameter  =  8  in..  Width  =  2  in.  -  Jones 
Beach  Sand  -  4°  Slope 
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Fig.  101  Run  #105  Transverse  Distribution  of  Normal  Stresses 

at  Central  Angles  16°  to  36°  Loose  Sand 
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Fig.  102  Run  #105  Transverse  Distribution  of  Normal  Stresses 

at  Central  Angles  -8*  to  12®  Loose  Sand 
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Fig.  103  Limiting  Normal  Stresses  due  to  Lateral  Failure 
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Fig.  104  Run  #107  Transverse  Distribution  of  Normal  Stresses 

at  Central  Angles  14*  to  22°  Dense  Sand 
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Fig.  105  Run  #107  Transverse  Distribution  of  Normal  Stresses 

at  Central  Angles  -8°  to  12°  Dense  Sand 


0  Distance  from  Center  of  Wheel  2in. 


Fig.  107  Run  #114  Transverse  Distribution  of  Normal  Stresses 

at  Central  Angles  -8°  to  12°  Dense  Sand 


Fig.  110  Orientation  of  Major  Principal  Stress  in  the  Front 
Field  Approximately  Corresponds  to  that  of  Instan¬ 
taneous  Velocity  Vector 
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Fig.  112  Pull  Coefficient  Versus  Predicted  Slip  Values  in 
Loose  Sand 
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01  770 
01  780 
01  790 
01  800 
01  81  0 
01  820 
01  830 
01  840 
01  850 
01  860 
01  87  0 
01  880 
01  890 
01  900 
01  91  0 
01  920 
01  930 
01  940 
01  950 
01  960 
01  970 
01  980 
01  990 
02000 
020  1  0 
02020 
02030 
02040 
02050 
02060 
02  0  7  0 
020  80 


T2=PI/2+. 5*(D4+DD1 )-AL 

QU1  =  (W0*COS ( EE1  >+CCl  )/COS (D2  ) 

QU2=QUP (D2,FF1  ) 

S IG=QU1 /QU2 

EP1  =EP0(FF1  ,T2,T1  ) 

S IG1 =S IG  *EP 1 
QU3  =  QUA CDD1 , FF1  ) 

QU4  =S I G 1 *QU3  *C  OS (DD1  )-CCl 
S IG2= (QJ  +  CC1  ) / (QU3  *C  OS (DD1  )) 

EP2  =  S  IG2/S  IG  1 
TF1  =TAN (FF1  ) 

P01  =ALOG  (EP2  ) 

THDEL=-P01  /  (2*TF1  ) 

AA1 =AL+THDEL 
IF  (QU4.GE.QJ)  GO  TO  340 

CALL  SLIP  (RR0,  EE1  ,  AA0  ,  GG0  ,  DD 1  ,  FF 1  ,CC0,DDF,  AL  ,  XI  ,  P2  ,  HH,  QQ 


JJ1  ) 

QD=Q J-OQ  (JJ1  ) 

IF  (OD.LT.0)  AX=AA0 
IF  (QD.GT.0)  AY=AA0 
IF  (AX*AY)  240  ,  250  , 240 

240  AA0= (AX+AY)/2 

IF  (ABS (QD)-.05*QJ)  310,310,245 
GO  TO  230 

245  GO  TO  230 

250  IF  (ABS (QD)-.05*QJ)  310,310,260 

260  IFCAA0.  EO.  ARE)  GO  TO  270 

GO  TO  290 

270  IF  (QQ  ( JJ1  ).  LT.  (I  -TOL)*S  IM)  GO  TO  280 

IF  (D0.  EQ.  DMAX)  GO  TO  810 
D0=B1 *D0  +QD*D0 /Q J+B4 
IF (D0. GT.DMAX)  D0=DMAX 
GO  TO  210 

280  S IM=QQ ( JJ 1 ) 

GO  TO  120 

290  AMD=AA0-AL 

QM=. 5* (Q J  +  OO  (JJ1  ) ) 

IF  (QD.LT.0)  ADD=-.02 

IF (QD*  GT • 0  )  ADD=.02 

AAD=AMD*QD/QM 

AA0=AA0+AAD+ADD 

IF(AA0. EQ. (AL+.01 ) )  GO  TO  310 

IF  (AA0.LT. AL)  AA0=AL+. 0 1 

IF  (AA0.GE.ARE)  GO  TO  300 

GO  TO  230 

300  AA0=ARE 

GO  TO  230 

310  DO  320  J=(M1-1  ),  (N1+1-JJ1  ),-l 

HK ( J )=HK ( J-N 1 + JJ1 ) 

QQ  (J  )=QQ ( J-Nl +JJ .  ) 

EE ( J )  =  EE (J-N 1 +JJ1  ) 
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020  90 

320 

CONT  INUE 

02100 

TORQ*0 

021  10 

LOAD=0 

02120 

DRAB**0 

02130 

DO  330  J=2*NN-1 ,N1 +1 -JJ1 ,-l 

02 1 40 

ARC=R0* (HH(J)-HH( J+l ) )/57.3 

02150 

QAV=. 5* (QQ (J+l ) +QQ (J) ) 

02160 

TAV=. 5*(EE(J+1  )  +EE  < J) ) 

02170 

AAV*=#  5 *  ( HH  (J  +  l  )  +HH  (J  )  )/57  •  3 

021  80 

LOA=  (TAV*S  IN  (AAV)+QAV*COS  (AAV)  )*ARC 

021  90 

DRA= (TAV*COS (AAV)-QAV*SIN (AAV) )*ARC 

02200 

TOR=R0*TAV*ARC 

022  1  0 

LOAD=LOAD+LOA 

02  2  2  0 

DRAB=DRAB+DRA 

02  2  3  0 

TORQ=TORQ+TOR 

02240 

330 

CONT  INUE 

02250 

, 

LOAD=B0*  (L0AD*C  OSD  (S0  )  +DRAB+S IND  (S0 

02260 

T  ORQB0  *T  ORQ 

02  2  7  0 

DRAB=B0* (DRAB-LOAD+S IND (S0 ) ) 

02  2  80 

GO  TO  410 

02  2  90 

340 

JJ1  =1  6 

02  3  00 

DALPH=57  #3*  (AAl  -AD/JJ1 

02310 

DO  350  J=(N1-1  ),  (N1+1-JJ1  ),-l 

02320 

HH ( J )  =  57 • 3  +AL+DALPH 

02  3  3  0 

QQ (J)=QJ 

02340 

EE ( J )  =  E J 

02350 

AL=HH ( J ) /57 • 3 

023  60 

IF  (J.  LT .  (N  1  -  1  )  )  GO  TO  350 

02  3  7  0 

3  50 

CONTINUE 

02  3  80 

AA0=HH(J)/57.3 

02  3  90 

T  ORQ=0 

02  4  0  0 

LOAD=0 

024  1  0 

DRAB=0 

024  2  0 

DO  360  J=2*NN-1 ,Nl +1 -JJ1 ,-l 

02430 

ARC=R0* (HH(J)-HH(J+1 ))/57.3 

02440 

QAV=.  5*  (QQ  (J  +  l  )  +QQ  ( J )  ) 

02450 

TAV=. 5* (EE (J+l  )  +  EE ( J ) ) 

02  4  60 

AAV=.5*(HH(J+1 )+HH(J) )/57.3 

02  4  7  0 

LOA=  (TAV*S  IN  (AAV )  +QAV+C OS  (AAV  ))*ARC 

024  80 

DRA= (TAV*C  OS (AAV ) -QAV+S IN (AAV ) ) *ARC 

024  90 

TOR=R0  *TAV*ARC 

02500 

LOAD=LOAD+LQA 

025  1  0 

DRAB = DRAB +DRA 

02520 

WE ( J )=B0* (LOAD*C OSD (S0 )+DRAB*S IND (S0 

02530 

DRB  (  J  )  =B0*  (DRAB-LOAD+S  IND  (S0  )  ) 

02540 

T  ORQ  =T  ORQ  +T  OR 

02  550 

TRQ (J)=B0*TORQ 

02  5  60 

IF  (J#GE#  (N  1  +1  )  )  GO  TO  360 

02  5  7  0 

IF  (WE(J).GT.L0>  GO  TO  370 

02  5  80 

IF  ( J.  LT •  (N  1  +1  )  )  GO  TO  370 
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02  590 

GO  TO  3  90 

02600 

3  60 

CONTINUE 

0261  0 

370 

WR= (WECJJ-L0 )/(WE(J)-WE(J+l  ) ) 

02  620 

LOAD=WE( J )-WR* (W  E  (  J  ) -WE ( J+  1  ) ) 

02630 

T 0RQ=TRQ ( J ) -WR* (TRQ (J)-TRQ (J  +  l  )) 

02  640 

DRAB=DRB(J)-WR*(DRB(J)-DRB(J+1  ) ) 

02  650 

AA0= (HHCJ)-WR* CHH(J)-KHCJ+1  )  ))/57. 

02  660 

I F  CA2  .  LE.  (5.73 *AA0  )  )  GO  TO  380 

02  670 

GO  TO  400 

02  680 

3  80 

A2=A2+B6 

02  690 

IF CA2. GT. AMAX)  A2=AMAX 

02700 

ITAL= 1 

02710 

GO  TO  410 

02720 

390 

LOAD=WE(J) 

02730 

DR AB=DRB ( J) 

02740 

TCRO=TRQ(J) 

02750 

GO  TC  410 

027  60 

400 

IFCD0. EQ.DMAX)  GO  TO  750 

02770 

KFR  =  1 

027  60 

410 

IF CK. EG.  CK3+K5) )  GOTO  830 

02  7  90 

IF  (K.  EO.  K3  >  GO  TO  940 

02  8  0  0 

420 

K=K+ 1 

02810 

LD (K)=LOAD 

02  82  0 

TR  CK ) =T CRQ 

02  830 

DR  CK)  =  DRAB 

02  840 

AR  CK)=A2 

02850 

DEL (K)=D0 

02  860 

AAE=57 . 3  *AA0 

02  870 

AE (K)=AAE 

02  880 

PU  (K ) =DR (K)/LD (K) 

02  890 

LC (K ) =LD CK)/L0 

02  90  0 

IF  (KFR.  EG.  1  )  GOTO  470 

02  91  0 

IFCITAL.EQ.  1  )  GO  TO  430 

02  920 

GO  TO  440 

02  930 

430 

IT  AL=0 

02  940 

GO  TO  210 

02  9  5  0 

440 

IF  CAR  (K).GE.  C0. 8+AECK)  )  )  GO  TO  450 

02960 

GO  TO  4  90 

02  970 

450 

IFCKA.EQ.  1  )  GOTO  460 

02  9  80 

AMAX=0. P+AECK) 

02  990 

KA=  1 

03000 

GO  TO  490 

03010 

4  60 

I F  C  AR  (K).GE.  (0.9*AECK)))  GO  TO  7  50 

03020 

GO  TO  490 

03030 

470 

IFCA2.EQ.AMIN)  GO  TO  4  80 

03040 

GO  TO  490 

030  50 

4  80 

IFCFU(K).  GE.  CPULC+TOP)  )  GO  TO  770 

030  60 

GO  TO  4  90 

03070 

490 

IF  CABS  CPU  CIO-PULC  >-TOP>  500,500,51 

030  80 

500 

I F CABS CLC CK ) - 1 >-TOL)  820,820,620 
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03090 

510 

IFCPUOO.GE.  (PULC+T0P>>  GO  TO  520 

03100 
031  10 
03120 

520 

IF  (PU(K).LE.  (PULC-TOP))  GO  TO  560 

GO  TO  620 

IF  (LC(K).GE.  l+TOP)  GO  TO  540 

03130 
03140 
031  50 
03160 
03170 

530 

IF  (LC(K).LE.  (1-T0P>>  GO  TO  550 
D0*B2*D0-B3* (PU (K)-PULC  5-B8 

IF(KFR.  EQ.  1  >  GOTO  530 

GO  TO  650 

A2=A2-B6 

031  80 
031  90 
03200 

540 

KFR=0 

GO  TO  650 

D0=B2*D0-B3* (PU (K)-PULC ) 

03210 

03220 

03230 

550 

A2=BR*A2-B8*LC  <K)-B8 

GO  TO  650 

D0=B2*D0-B3*  (PU  (JO-PULC  >-B8 

03240 

03  2  5  0 
03260 

560 

A2»B1*A2+B5*LC (K) +B4 

GO  TO  650 

IF(KFR.EQ. 1 )  GO  TO  570 

03  2  7  0 

03  2  80 

03  2  90 
03300 
03310 
03320 

570 

IF  (LC(K).GT.(1+T0L>>  GO  TO  600 
IF(LC(K).LT.  (1-T0L)>  GO  TO  610 

D0=B1 *D0-B3* (PU (K)-PULC ) +B4 

IF(KFR.EO.l)  GO  TO  570 

GO  TO  650 

IF(A2.GE.  (AMAX-1  )  )  GO  TO  750 

03330 

03340 

03350 

580 

IF(A2.GE.43*AA0>  G  0  TO  580 

GO  TO  590 

D0=B1 *D0-»  5*B3* (PU (K)-PULC )  +  • 5*BT 

03  3  60 

590 

A2=A2+B6 

03370 
033  80 

600 

GO  TO  650 

D0=B1 *D0-B3* (PU (K) -PULC ) +B4 

03  3  90 
03400 
0341  0 

610 

A2*B2*A2-0. 5*B5*LC (KJ-B8 

GO  TO  650 

D0=B1 *D0-B3* (PU (K)-PULC >+B4 

03  4  2  0 

03  4  3  0 

03  4  4  0 

620 

A2=B1 *A2+B5*LC (K)+B4 

GO  TO  650 

IF (LC (K).GT.  1 >  GO  TO  630 

03  4  50 

03  4  60 

630 

IF(LC  (K).LT.  1  )G0  TO  640 

A2  =B2  *A2 -B  5*LC (K) 

03  4  7  0 
034  80 

03  4  9(3 

640 

I F  (A2  •  LT  .  AM  IN  )  A2=AMIN 

GO  TO  650 

A2=B1 +A2+B5/LC (K) 

03  5  0  0 
03510 

03  520 

650 

I  F  (A2 . GT . AMAX >  A2=AMAX 

GO  TO  650 

IF  (A2.LT.AMIN)  A2=AM IN 

03530 

03  540 
03550 

03  5  60 

03  5  7  0 

03  580 

660 

IF (A2»  GT .AMAX)  A2=AMAX 

IF  (D0 • LT . DM IN  )  D0=DMIN 

IF  (D0.GT.DMAX)  D0=DMAX 

I F ( KFR •  EQ •  1  )  GOTO  660 

GO  TO  670 

KFR=0 
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03  5  90 

GO  TO  210 

03  600 

670 

IF (AR (K-l  ). EQ»  AR  CK) )  GO  TO 

680 

036  1  0 

GO  TO  800 

03  6  2  0 

680 

IF  (PU (K) • GE.PULC )  GO  TO  690 

03  630 

GO  TO  700 

03  640 

690 

IF  (PU(K~1 ).LE.PULC)  GOTO 

710 

03  650 

GO  TO  720 

03  660 

700 

IFCPUCK-l  ). GE.PULC  )  GO  TO  ‘ 

710 

03670 

GO  TO  720 

03  680 

710 

D0  = (DEL (K>+DEL(K-1  ))/2 

03690 

IF  (D0.LT.DMIN)  D0=DM IN 

03700 

IF  (D0.GT.DMAX)  D0=DMAX 

03710 

GO  TO  800 

03  7  2  0 

720 

IF  (DEL (K-l ).EQ»DEL(K))  GO 

TO  730 

03730 

GO  TO  800 

03740 

730 

IF(A2. EQ.AMAX)  GO  TO  740 

03750 

IF(A2.  EQ.  AMIN  )  GO  TO  7  60 

037  60 

IF(D0.  EQ.  DMAX)  GO  TO  7  80 

03770 

IF  (D0.EQ.DMIN)  GO  TO  7  9 0 

037  80 

GO  TO  800 

037  90 

GO  TO  820 

03  800 

740 

IF  (LC  (K).LE.  (1 -TOL)  )  GO  TO  750 

03  81  0 

GO  TO  820 

03  820 

7  50 

SLP= 1 .0 

03  830 

GO  TO  960 

03  840 

7  60 

IF  (LC  (K).  GT.  (1 +TOL)  )  GO  TO 

770 

03  850 

GO  TO  820 

03  860 

770 

SN  K=0 

03  870 

GO  TO  960 

03  880 

GO  TO  820 

03  890 

7  80 

IF (PU (K).LE. (PULC-TOP ) )  GO 

TO  750 

03  90  0 

GO  TO  820 

03  91  0 

GO  TO  820 

03  9  2  0 

7  90 

IF(PU(K).GE.  (PULC+TOP)  )  GO 

TO  770 

03  930 

GO  TO  820 

03  9  4  0 

GO  TO  820 

03  9  50 

800 

GO  TO  210 

03960 

810 

F5=F2/57 . 3 

03  970 

TA2=TAN (F5 ) 

03  980 

PS  I=C2/TA2 

03  990 

F2  =B7  *F2 

04  0  0  0 

F5=F2/57.3 

0401  0 

TA2=TAN (F5 ) 

04020 

C2  =PS I *TA2 

04030 

IF  (KOC  *  G  E.  3  )  SNK=0 

04  0  4  0 

GO  TO  960 

04  0  5  0 

KOC  =  KOC  +  l 

04  0  6  0 

GO  TO  190 

04070 

820 

GO  TO  942 

040  80 

830 

K6=K3+K5 
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04090 
04100 
04110 
04120 
04130 
04140 
041  50 
041  60 
04170 
041  80 
041  90 
042  0  0 
042  1  0 
04  2  2  0 
04  2  3  0 
04  2  4  0 
04  2  50 
04260 
04  2  7  0 
04  2  80 
04290 
04  3  0  0 
04310 
04  3  2  0 
04  3  3  0 
04340 
043  50 
04  3  60 
04370 
043  80 
04  3  90 
04  4  0  0 
044  1  0 
044  2  0 
04  4  3  0 
04  4  4  0 
044  50 
04460 
044  7  0 
04480 
044  90 
04500 
0451  0 
04520 
04530 
04  540 
04  5  50 
04  5  60 
04570 
04  5  80 


LKWH.SRC  09- JAN -7 3  08*57 

DO  840  K= 1 » K6 
IF(LC  (K).GT.  1  .0)  GO  TO  850 
840  CONTINUE 

GO  TO  900 
850  K1=K 

LMX=LC  <K1 ) 

860  DO  870  K=1,K6 

IF  (LC  (K) • LT* 1 « 0  )  GO  TO  880 
870  CONTINUE 

DO  87  5  K=1,K6 
IF(LC CK). EQ. 1 .0)  GO  TO  885 
87  5  CONTINUE 

880  K2  =  K 

LMN=LC (K2 ) 

GO  TO  905 
885  K4=K 

IF(K1.EQ.0>  GO  TO  890 
IF(K2#EQ»0)  GO  TO  895 
GO  TO  905 
890  LMX=LC(K4) 

K1  «=K4 
GO  TO  90  5 
895  LMN=LC(K4) 

K2  =  K4 
GO  TO  905 
900  GO  TO  750 

GO  TO  940 

90  5  DO  915  K=1<K6 

IF(LC  (K).LT.LMX)  GO  TO  910 
GO  TO  915 
910  LMX=LC (K) 

915  CONTINUE 

DO  920  K= 1 t K6 
IFCLC (K).EQ.LMX)  GO  TO  925 
920  CONTINUE 

925  KM=K 

DO  930  K= 1 j K6 
IF  (LC (K).GT.LMN)  GO  TO  98 
GO  TO  930 
928  LMN  =LC (K) 

930  CONTINUE 

DO  932  K- 1 j K6 

IF  (LC  (K).EQ.LMN)  GO  TO  934 
932  CONTINUE 

934  KN  =  K 

IF  (PU  (KN  )  •  LE*  PULC  )  GO  TO  936 
GO  TO  980 

936  I F  (PU  (KM  )  •  G  E«  PULC  )  GO  TO  938 

GO  TO  986 

93  8  PUR= (PU (KN )-PULC  >/ (PU (KN  )-PU (KM) ) 
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04590 
04600 
04610 
04620 
04630 
04640 
04650 
04660 
0T670 
04680 
04690 
04700 
047  1  0 
04720 
047  3  0 
04740 
047  50 
047  60 
047  7  0 
04  7  80 
047  90 
04  800 
04810 
04  820 
04830 
04  840 
04  850 
04  860 
04870 
04  8  80 
04  890 
04900 
04910 
04920 
04  930 
04  94  0 
04  950 

LSL) 
04  9  60 
04  970 
04  9  80 
04990 
05000 
05010 
05020 
05030 
0504  0 
05050 
050  60 
05070 
050  80 
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K=K6+1 

LD (K ) =LD (KN  }  +PUR* <LD (KM ) -LD (KN ) ) 

TR  (K)=TR (KN ) +PUR* (TR (KM )-TR (KN ) ) 

DR(K)=DR(KN  )+PUR*(DR(KM)-DR(KN)  ) 

AAE«AE(KN )+PUR* (AE(KM)-AE(KN ) ) 

A2=AR (KN  J+PUR* (AR (KM)-AR(KN ) ) 

D0=DEL (KN )+PUR* (DEL(KM)-DEL (KN ) ) 

LC (K)*LD(K)/L0 

IF(LC(K).GE.  (1+TOP))  GO  TO  984 
IF(LC(K).LE.  (1-TOP))  GO  TO  980 
GO  TO  820 
940  B1  =  1 • 03 

B3  =  1  5 
B4=0.5 
B5=  1 
B6=  1 

GO  TO  420 

942  SNK=R0* ( 1 -COS (AA0 ) )*1 2 

SP=  (1  -D0/DMAX) 

IF  (SP.  EQ « 0  )  GO  TO  944 
SLP=SLJ-SLK*ALOG (SP) 

GO  TO  960 
944  SLP= 1 • 0 

GO  TO  980 
950  FI =F0/57 • 3 

TA3=TAN (FI  ) 

C0=C0+0. 5*SIM*TA3 

RR*R0 

BB=B0 

LA=L0 

DB0=DB 

SS=S0 

CC=C0 

GG=G0 

ZJ=SLJ 

ZK=SLK 

CALL  LKWC  (RR,BB,LA,DB0,SS,CC,GG,ZJ,ZK,TOL,TOP,LL,DP,TQ,S 


K=  1 

LD (K)=LL 
DR (K)=DP 
TR (K)=TQ 
SNK=S  I 
SLP=SL 
GO  TO  960 

960  IF(SLP.EQ.  1 .0)  GOTO  980 

IF  (SNK.  EQ»  0  )  GO  TO  984 
PRINT  965, LD (K)  ,  DR (K) 

965  FORMAT  (1H  , ’LOAD*  ' ,F10.2,’  DRAWBAR  PULL  =  ',F10.2) 

PRINT  970,  TR  (K) 

970  FORMAT  (1H  ,  *TORQUE=  • ,F10.2) 
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050  90 

PRINT  97  5,SNK,SLP 

05100 

97  5 

FORMAT  C1H  ,’SINKAGE=  *,F10.2,*  IN  SLIP*  »,F10.3> 

05110 

GO  TO  999 

05120 

980 

PRINT  982 

05130 

982 

FORMAT  C 1 H  »'  NO  GO  :  100  X  SLIP  ’  ) 

05140 

GO  TO  999 

051  50 

984 

PRINT  986 

051  60 

986 

FORMAT  (1H  » ’USE  HARD  SURFACE  FORMULA  OR  VARIABLE  DELTA* 

) 

05170 

990 

CALL  RTIME  (11 ) 

051  80 

PRINT  994/  11 

05190 

994 

FORMAT  (1H0/*TIME*  * ,110) 

05200 

999 

STOP 

SYSTEM?.  . 
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00100 

SUBROUTINE  SLIP  <R0,E1  *  A0*  G0>  D1  ,  FI  *  C0..  DF,  AM,  XT,  PL/ P,  Q  ,  E,  J 

1  > 

00110 

DIMENS  ION  X(4  9,  17  ),  Z  (49,  1  7  >,S  <49,  17>,T(49,  17  ) 

00120 

DIMENSION  A(35),B(35>,C<35> 

00130 

DIMENSION  D  (35  ) 

00  1  40 

DIMENSION  H(35),Q  (35),E(35> 

001  50 

IMPLICIT  REAL  (L) 

001  60 

DEL(D9,F9>=ASIN (SIN (D9)/S IN (F9 )  ) 

00170 

QUA(D9,F9)=C0S (D9J+SQRT (COS (D9>**2-C0S (F9>**2  ) 

001  80 

QUP(D9,F9)=C0S (D9 ) -SORT (COS (D9>**2 -COS (F9)**2) 

001  90 

EP0(F9,T9,T8)  =  EXP  (2*  (T9-T8>*SIN  (F9)/C0S  (F9  ) ) 

002  0  0 

TAN(F9)=SIN(F9)/C0S (F9) 

002  1  0 

DIS  (A9,D9,T9,T8,F9)  =  R0*(A9-D9)*EXP(  (T9-T8)*SIN  (F9)/C0S  (F 

9)  ) 

0022  0 

M=4  9 

00230 

N=  1  7 

00240 

PI =3. 1 41 59 

00250 

1  6 

W0  =  0 

002  60 

TF=TAN (FI ) 

00  2  7  0 

DAF=TAN (D1  ) 

00  2  80 

11  =2*N-l 

00290 

C 1 =C0/TF 

00  3  0  0 

U1  =P  I  /4  -  FI  /2 

00310 

F3  =  1  -S  IN  (FI  > 

00  3  2  0 

V3  =S  IN  (El  -FI  >/COS  (FI  ) 

00330 

V4=C  CS (El -FI  )/COS (FI  ) 

00340 

V5=S IN (El +F1 )/COS (FI ) 

00  3  50 

V6=C OS  (El  +F1  )/CCS  (FI  ) 

00360 

D2=ATAN  (W0*S  IN  (El  )/ (Cl  4W0*COS  (FI  )  )  ) 

00  3  7  0 

D3=DEL(D2,F1  ) 

003  80 

Tl=(D3-D2)/2 

00390 

D5=D1  *(1  4DF5/2 

00400 

D4  =  DEL(D5,F1  ) 

00410 

T2=PI/24. 5* (D4+D1 ) -A0 

00420 

W1  =  (W0*COS(E1  )  4  C  1  >/C0S(D2) 

00  4  3  0 

W2=QUP(D2,F1 ) 

00440 

S I =W 1 /W2 

004  50 

T3=. 75*P 14. 5*F1 -T24T1 -A0 

00460 

LI =DIS (A0,AM,T2,T1 ,F1 ) 

00  4  70 

L2  =2 . 5  *L 1 *COS (U1  )*COS (T3-F1 )/C OS (FI  ) 

004  80 

L=L2 

00  4  90 

20 

DO  40  J=1  >N 

00  5  0  0 

I=N4l -J 

0051  0 

A  J  =  J 

00  520 

AN=N 

00530 

Z ( I, J>=0 

00540 

S (I, J)=51 

00  5  5  0 

T (  lj J )=T 1 

00  5  60 

X( I, J>= (AJ-2 )  *L/ (AN -2 ) 

00  5  7  0 

X(N,  1)=0 

00580 

X(N-1,2)=. 5*L/ (AN -2  ) 

00  5  90 

40 

CONTINUE 

00600 

50 

DO  70  I  =N j  (2*N-1  ) 

0061  0 

J=1 
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00  6  3  0 

A I “FLOAT (I) 

00640 

AN  =  FL  OAT  (N  ) 

00650 

XCI,J>=0 

00  660 

Z(Jj  J)=0 

00  670 

T  < I *  J  )=T 1 + (T2-T 1  )* (A I -AN )/ (AN-1  ) 

00660 

POV=EPO(Fl,TU,J),Tl  ) 

00  6  90 

S(I,J>=POV*S  (N,l  ) 

007  0  0 

70 

CONTINUE 

00710 

Q2=QUA  CD1  ,F1  ) 

00  7  2  0 

Q  (1  )=Q2*S  (2*N-1  ,  1  >*COS  CD1  >-Cl 

00  7  3  0 

EC1  )=  (Q  (1  )+Cl  ) *DAF 

007  40 

H  ( 1  )  =  57.3*A0 

00750 

DO  200  J=2,N 

00  7  60 

80 

DO  180  1=  (N+2-J),  (J+2*(N-1  >) 

00  7  7  0 

IF  Cl.  EQ.  <J+2*(N-1  )  )  )  GO  TO  130 

00  7  80 

K=0 

007  90 

TH1=TCI, J-l )4U1 

00  800 

TH2=T (1-1 ,J>-U1 

0081  0 

SI1=S  (I.J-1  ) 

00  820 

SI2=S(I-1*J) 

00  83  0 

V7  =2  *S  C  I  - 1  ,  J  )  *S  (  I ,  J-  1  ) 

00  840 

U9=S ( I - 1  , J>  +  S (I, J-l  ) 

00  8  50 

V8=(T(I,J-1  >-T(I-l,J)>*TF 

00  860 

V9=2*S  (1-1 , J)*S  Cl, J-l )*V8 

00  870 

U6=2*TF*<S (I, J-l  >*T  <1, J-l )+S (I- l, J)*T  <1-1 , J) > 

00  880 

90 

VI  =TAN  <thi  > 

00  890 

V2=TAN  CTH2  ) 

00  900 

XI=V1 *XCI, J-l ) 

0091  0 

XJ=V2*X< 1-1  ,  J) 

00  92  0 

VI 2=1 /(VI -V2 ) 

00  93  0 

X( I, J)=V12* (Z ( 1-1 , J)-Z (I, J-l  )+XI-XJ> 

00940 

Z( I, J)=Z(I-1 . J)  +  CXCI, J)-X(I-1  ,  J) )*V2 

00  950 

AA=V3* (X(I, J  >-X(I-l >  J ) >+V4*(Z(I, J)-Z (1-1 , J) } 

00  960 

BB=V5*(X(I>J)-X(I> J-l ) >+V6*(Z(I, J)-Z (1, J-l )> 

00  9  7  0 

U5=S (  Ij J-l  ) - S ( 1-1 , J) 

00  980 

IF  (U9.LE.0)  GO  TO  110 

00  9  90 

GO  TO  120 

01  000 

1  10 

PRINT  112 

01010 

1  12 

FORMAT  C1H-*' CANT  COMPUTE  CASE  s  CHECK  INPUT •  > 

01040 

120 

S  (I, J)=CV7+V9+G0*(SI1*AA+SI2*BB> )/U9 

010  50  ■ 

T (I, J>  =  CU5+U6+G0*(BB-AA) >/(2*TF*U9) 

010  60 

IF  (  K.  EQ. 1  )  GO  TO  1  80 

01  070 

TH1=.5*(T(I,J-1  )  +T  (  I .  J  >  )  +U1 

010  80 

TH2  =  .5*(T(I-1,J>+T(I,J))-U1 

01090 

S 1 1  =.  5*  (S  <  1/  J-  1  >+S  (I,  J)  ) 

01  100 

SI2  =  .  5*  (S  <1-1  ,  J)+S  Cl#  J)  ) 

01  1  10 

U6=2*TF*(SI1  *T  (I, J-l  )+S  I2*T  C  1-1  ,  J)  ) 

01  120 

V7=S  11  *S  (1-1  ,  J)+S  I2*S  (I,  J-l  ) 

01  130 

V9=2*S  1 1  *SI2*V8 

01  140 

U9=S  1 1 +S  12 

01  1  50 

K=  1 

01  1  60 

GO  TO  90 
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01  170 

130 

CH=R0  *S IN (A0 > 

01  1  80 

CHl=CH+y(J+2*N-3.»J-l  ) 

01  1  90 

A  (  J )=AS IN (CHI /R0 ) 

01200 

D  (  J  )  =D  1  *  ( 1  -  (1  -DF>*  (A  (J)-A  (2)  )/  (AM-A(2  )  )  ) 

01  21  0 

DA J  =  TAN (D  ( J )  ) 

01220 

B1=DEL(D(J),F1  > 

01230 

B2  =  . 5*  CD ( J )+Bl  ) 

01  240 

T  (Ij  1-2*  (N-l  )  )  =P  I/2+B2-ACJ) 

01  250 

TH3  =  .  5*  (T  (  I,  1-2*  (N-l  >  >+T  (  I  -  1  ,  I -2  *  (N- 1  )  )  ) 

012  60 

TH4  =T  H3 -U 1 

01270 

TH5  =  TAN  (TH4  ) 

01  280 

TAJ=TAN  (A  ( J)  ) 

01  290 

Z1  =1  /TH5 

01  300 

Z2  =  Z  1 *TAJ 

01  3  1  0 

Z3  =1/(1  +  Z2  ) 

01  320 

Z4  =  Z1  *  (7  (1-1  ,  J-l  >-Z ( 1-1 , J  )  > 

01  330 

X( Ij J)=Z3*(X(I-1, J>+Z2*X(I-1 , J-l )+Z4 ) 

01  340 

Z(I*J>=Z(I-1,J-1 >*TAJ* (X(I-1 ,J-1 )-X( I, J> ) 

01350 

A A  =  V3*(X(I>J)-X(I-1 , J> )+V4*(Z (I, J)-Z ( 1-1 , J) > 

01  360 

U3=2*S ( 1-1 , J)*TF* (T ( Ij J)-T ( 1-1 , J ) ) 

01  370 

S  (  I*  J  )=S  (  I  -  1 J)+U3+G0*AA 

01  3  80 

Q  1  =  QUA (D  ( J  )  >  FI  ) 

01  390 

Q  ( J  )=Q  1  *S  ( I  ,  J  >*C  OS  (D  (J  )  )  -C  1 

01400 

E  ( J  )= (Q (J)+C1  ) *DA J 

01410 

CH2=CH+X( 1, 1-2* (N-l ) ) 

01420 

C H3  =AS  IN  (CH2/R0  ) 

01  430 

F  ( J  )  =  57 . 3  *C  H3 

01  440 

IF  (XX.  EQ.  1  )  GOTO  155 

01450 

AN=DF*ATAN ( (Q ( J  ) +PL  >  *DA  J/Q (J)  ) 

01460 

155 

IF  (CH3  .GE.  (AM+.002  )  )  GO  TO  176 

01470 

IF  (CK3.  LE.  (AM-.  002  )  )  GO  TO  160 

01480 

IF  ( J.  LE.  7  )  GO  TO  1  65 

01490 

GO  TO  220 

01  500 

1  60 

I F  (  J.GE.7)  GO  TO  170 

01  510 

1  65 

L= • 5*L+. 05 

01  520 

GO  TO  20 

01  530 

1  70 

HI  =K(J)-H(J-1  ) 

01  540 

H2  =  57 .3 *AK-K (J-l  ) 

01  550 

H3=H2/H1 

01  560 

H4  =  H3*(X(N  +  1  -J,  J)-X(N-J  +  2,  J-l  )  ) 

01  570 

X(N  +  1  -J,  J)=X(N-J+2,  J-l  )  +H4 

01  580 

GO  TO  80 

01  590 

176 

IF  (J.LT.N  )  GO  TO  1  80 

01  600 

L= 1 . 5*L+. 1 

01  61  0 

GO  TO  20 

01  620 

1  80 

CONT  INUE 

01  630 

200 

CONT  INUE 

01  640 

220 

J1  =J 

01  650 

RETURN 

01  660 

1  000 

END 
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00  100 

' 

SUBROUT  INE  LKWC  (R0  » B0  »  L0*  DB*  S0  *C0j  G0  j  S  J.*  SKjTL » TP*  LL*  DP>T 

Q,SI,SL) 

00  1  1  0 

DIMENSION  HH<35)/QQ(35)>EE(35) 

00120 

DIMENSION  WE  <3  5  )#DRB  (3  5)>TRQ  (3  5) 

001  30 

D IMENS  ION  LD  ( 1  5  )  ,  DR  (  1  5  ) , TR  <  1  5 )  #AR  <  1  5  ) * DEL (15) 

00140 

D I MENS  I  ON  PU  ( 1  5  )  *  LC  ( 1  5  ) 

001  50 

DIMENSION  AE (15) 

001  60 

IMPLICIT  REAL  (L) 

00170 

TAN CF9  )=S IN (F9)/C0S CF9) 

001  80 

P 1  =  3  » 14159 

001  90 

DEAMAX= 1 . 57 

00200 

DEAM  IN  = .  02 

00210 

DB 1 =DB+L0*S IND (S0 ) 

00220 

PULC=DB1 /L0 

00  2  3  0 

DEA=AS IN C2*PULC  +  • 1  ) 

00240 

IF  CDEA.GT.DEAMAX)  DEA=DEAMAX 

00  2  5  0 

IF(DEA.LT.DEAMIN)  DEA=DEAM IN 

00260 

10 

AMAX=21 

00270 

AM  IN  =  1 

00  2  80 

ARMIN=-AMIN/57 .3  . 

00  2  90 

20 

ALF-DEA/2 

00  3  0  0 

A2=-ALF/3 

00310 

30 

DO  60  J=l,35 

00  320 

IF  CJ.GT. 11)  GO  TO  40 

00330 

HHCJ)=A2+ CJ-1 )*4*ALF/30 

00  3  4  0 

QQ  CJ)=C0*  Cl  +COS  CDEA )  +P I -ALF+HH  ( J  )  ) 

003  50 

EE(J)  =  C0*S  IN  CDEA) 

00360 

IF  CJ.LE.  1  1  )  GO  TO  50 

00  3  7  0 

40 

HH  C J )=ALF ♦ CJ-1 1  ) *2/57. 3 

00380 

QQ  CJ)=C0*Cl  +COS  CDEA)4PI+ALF-HHCJ)  ) 

00  3  90 

EE  C  J  )=C0*S  IN  CDEA) 

00  4  0  0 

50 

AL=57.3*HH<J) 

00410 

60 

CONTINUE 

00  4  2  0 

70 

TORQ=0 

00430 

LOAD=0 

00440 

DRAB=0 

00450 

DO  65  J=2,35 

004  60 

ARC=R0*CHH(J)-HHCJ-1 )) 

00  4  7  0 

QAV=. 5*  CQQ  CJ-1 )+QQ  CJ) ) 

004  80 

TAV=. 5* CEE CJ-1 )+EECJ)) 

00490 

AAV=. 5* CHHCJ-1  )+HHCJ)) 

00  5  0  0 

LOA=  CTAV+S  IN  CAAV  )  +QAV*C  OS  CAAV)  )*ARC 

0051  0 

DRA=  CTAV+COS  CAAV ) -QAV*S IN  CAAV)  )*ARC 

00  520 

TOR=R0*TAV*ARC 

00  530 

LOAD=LOAD+LOA 

00  5  4  0 

DRAB=DRAB  +DRA 

00550 

WE  C  J  )  =B0  *  CLOAD+C  OSD  CS0  )+DRAB*S  IND  CS0  )  ) 

00560 

DRB C J )=B0* CDRAB-LOAD*S IND CS0  )  ) 

00  57  0 

TORQ=TORQ+TOR 

00580 

TRQ C J )=B0*T  ORQ 

00590 

IF  CWECJ).GT.L 0)  GO  TO  80 

00  600 

65 

CONTINUE 
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00610 

IF  (WE(J).LT.  (1-TL)*L0)  GO  TC  150 

00  62  0 

GO  TC  80 

00  630 

80 

I F  CHK ( J )  .  LT .  (-2  *A2  )  )  GO  TO  90 

00  640 

GO  TO  100 

00  650 

90 

IF  (A2.GE.  ARMIN  )  GO  TO  130 

00  660 

A2=A2  + . 02 

00  6  7  0 

IF  (A2.  GE.  ARMIN  )  A2=ARM  IN 

00  6  80 

GO  TC  30 

00  690 

1  00 

WR=  (WE(J)-L0>/(WF.(J)-WECJ-1  )  ) 

00  7  0  0 

LOAD  =  WE  ( J  )  -VR*  (WE(J)-WE(J-l  )  > 

00710 

T ORG  =TRQ  (J)-VR*  CTRQ  (J)-TRG  C J-l  )  ) 

00720 

DRAB=DR8  C  J ) - VR* (DRB (J)-DRB C  J  -  1  ) ) 

00  7  3  0 

AA0  =  CHHC  J  )-WR*  CHHCJ  )-HHCJ- 1  )  ;  ) 

00740 

1  10 

IFCK.EQ.  1  5  >  GO  TO  200 

00  7  5  0 

K=K+  1 

007  60 

LD  CK)=LOAD 

00  7  7  0 

TR  (K)=TGRQ 

007  80 

DR CK )  =  DRAB 

00  7  90 

AR  CK )=A2 

00  800 

ARK=57 . 5 *A2 

00  81  0 

DEL  CK)  =  DEA 

00  820 

AA  E=  57 . 3 *AA0 

00  830 

AE  CK)=AAE 

00  840 

PU (K)=DR CK ) /LD (K) 

00  850 

LC  ( K )  =  1.D  CK)/L0 

00  860 

IFCA3S  C*V  (K)-PULC  )-TP)  120 . 120, 1  50 

00  870 

1  20 

LL=LD (K ) 

00880 

DP=DR (K) 

00  890 

TQ=TR  GG 

00  9  0  0 

S  I  =R0*  (1  -COS  (AA0  )  )*1  2 

00  91  0 

SP  =  C 1 -DEA/DFAMAy ) 

00  920 

IF  <SP.  E9.0  >  GO  TO  140 

00  93  0 

SL=S J-SK*ALOG (SP) 

00  940 

GC  TC  200 

00  950 

130 

S  1  =  0 

i?0  960 

SL  =  0 

00970 

GC  TO  200 

itfi  9  80 

1  40 

SL= 1 .0 

'0  990 

GC  TC  200 

01  000 

1  50 

IF  CPU  CTO  .GE-  ''f  JLC+TP))  GC  TC  •  CO¬ 

01  01  0 

IF  f Pi  K’i.LL.  CPULC-TP)  )  GO  TO  170 

01020 

1  60 

IFCDEA.EO.DHAMIM  )  GC  TO  130 

01  030 

DEA  =  DEA  +  CPUl.C-PU  CK)  )-0. 0  5 

01  040 

I  FCDEA.LT.  DEAN  TN  >  DEA=DFAM  IN 

01  0  50 

GO  TO  1* 

01  060 

1  70 

IF  (DEA.  EQ-DEAMAV)  GO  TO  140 

01  070 

PEA  =  DFA*  CPT,1  C-PU  CK)  )*0.  n-'J 

01  0  80 

IFCDEA-GT.  DEAMA  X )  D  EA  -  DEAN  A  X 

01  090 

go  ~o  in 

01  100 

200 

RETURN 

01  1  1  0 

STOP 
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